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ABSTRACT 


EVALUATION OF EXERCISE-RESPIRATORY 
SYSTEM MODIFICATIONS AND INTEGRATION 
SCHEMES FOR PHYSIOLOGICAL SYSTEMS 

R.R. GALLAGHER, Ph.D. 

DEPARTMENT OF ELECTRICAL ENGINEERING 
KANSAS STATE UNIVERSITY 
MANHATTAN, KANSAS 

Exercise subroutine modifications are implemented in an exer- 
cise-respiratory system model yielding improvement of system res- 
ponse to exercise forcings. A more physi ol ogi cally desirable res- 
piratory ventilation rate in addition to an improved regulation of 
arterial gas tensions and cerebral blood flow is observed. A res- 
piratory frequency expression is proposed which would be appropriate 
as an interfacing element of the respiratory-pulsatile cardiovas- 
cular system. 

Presentation of a circulatory-respiratory system integration 
scheme along with its computer program listing is given. The inte- 
grated system responds to exercise stimulation for both nonstressed 
and stressed physiological states. Other integration possibilities 
are discussed with respect to the respiratory, pulsatile cardiovas- 
cular, thermoregu! atory , and the long-term -circulatory -systems . 
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1. INTRODUCTION 
1 . 1 Overv i ew of Model i ng Effort 

The major objectives and emphasis of the overall physiological 
system modeling effort focus on the development of an effective 
whole-body algorithm as being the most significant feature. 

(1) Utilizing system theory, mathematical descriptions of physio- 
logical functions, appropriate empirical interrelationships bet- 
ween variables, and a computer implementation and display system 
complex physiological systems have been realized. 

Physiological systems that were deemed important to the pro- 
ject's overall goals included the respiratory, fluid-electrolyte 
balance circulatory, pulsatile cardiovascular, and thermoregulatory 
systems. Proper evaluation, modification, and adaptation of 
selected models (2-5) have been made with their implementation co- 
ordinated with the project's goals. The indicated references pro- 
vide the necessary details of the models. Therefore, only specific 
modifications will be discussed in this report. 

The computer program for the integrated circulatory and res- 
piratory control systems which will be presented in more detail 
in later sections is included in the Appendix. This program is 
comment laden and provides an adequate description of the involve- 
ment of each subroutine within itself and also its relationship 
to other subroutines in the program. Many of the comments in the 
respiratory system component refer to Grodins' basic model. (2) 
Block diagrams and explanatory material for each of the two major 
systems (circulatory and respiratory) are given by White (6) and 
Gallagher (7). 
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Since the primary thrust of the modeling revolves around the 
effects on parameters that are influenced by the space environment, 
variations in gravitational forces, external body pressures, inspired 
gas composition, and thermal conditions are some of the inputs 
utilized in determining dynamic physiological responses and asso- 
ciated hypotheses. Implementation and coordination of the simulations 
with major Skylab experiments provide the research-application im- 
petus. Such a research effort demands a well-defined set of tasks 
which are fulfilled in sequence and/or concert. Some of the major 
milestones which are detailed by numerous- sped fi c tasks include 

(a) evaluation and modification of major physiological sub- 
systems , 

(b) development and evaluation of hypotheses regarding physio- 
logical systems responses' to stresses and altered en- 

vi ronments , 

(c) implementation of integrated system simulations and verifi- 
cation of responses, and 

(d) implementation of supporting software systems for col- 
lecting, storing, and correlating simulated a nd e xperi - 
mental data. 

With utilization of this complete whole-body simulation or its 
selected subsystems, revised techniques for remote medical care 
can be evaluated and implementation procedures established for 
future manned spaceflights. In addition, potential diagnostic 
and therapeutic capabilities may emerge from correlation of simula- 
tion and physiological testing during real time medical monitoring. 

In keeping with the spirit of the project one should not lose 
sight of the important contribution that simulation provides in 
understanding physiological processes under the influence of normal 
and altered environments. 

1 .2 Exercise and the Respiratory System Model 

The respiratory control system is an excellent example of a 
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physi ol ogical system that illustrates the influence of exercise 
upon its variables' responses. The magnitude and length of the 
exercise stimulus as well as the physical condition of the sub- 
ject play integral roles in determination of the system's responses. 
One aspect of exercise simulation discussed in the following report 
involves the .respiratory control system functioning as an indepen- 
dent system. Modifications of the exercise subroutine are eval- 
uated illustrating improved control of ventilation rates and arte 
i a 1 and compa rtmental gas tensions. 

1 . 3 Integration of Physi ol ogi ca,l System Models 

There are several approaches which can be taken in the process 
of integrating the simulations of the respiratory and circulatory 
systems. Since both systems' dynamic responses are influenced by 
exercise and with exercise variations and their associated physio- 
logical forcing roles being major components of the manned space- 
flight medical experiments, the integration scheme revojves around 
the implementation of exercise simulations. As implied in Section 
3 there are extensions in the utility of these models for non- 
spaceflight experiments. (8) 

In the following report emphasis is placed upon the inter-, 
face component of the integrated model. For details concerning 
other variable formulations in the individual modes one should re- 
fer to the references. (2, 3, 6, 7). The interface allows for a 
minimal amount of alteration of each existing system while im- 
proving or supplementing the simulation responses of each system. 

The integrated system can be described in terms of three 
functional components; stimuli, basic control systems, and the 
interfacing system as illustrated in Figure 1. Major categories 
of stimuli include exercise, parameter, and environmental variations. 




SXIMUL1 


BASIC CONTROL SYSTEMS 


INTERFACING SYSTEM 



Figure 1 


Integration scheme for respiratory and circulatory 
control systems. 
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As previously mentioned emphasis is placed upon the exercise stim- 
ulus since several manned spaceflight experiments stress 'man's 
physiological performance under various exercise levels in an al- 
tered environment. 

Since all of the basic control system's parameters are not 
common to both subsystems another stimulus includes the option of 
changing subsystem parameters and properties. This allows for in- 
vestigation of the response of one system to a particular physio- 
logical change when that response is only indirectly related to that 
change. Therefore, another dimension of complexity has been added 
to physiological system modeling which has been previously un- 
a vai Table. 

Another stimulus which directly affects the respiratory system 
model is the gaseous composition of the environment. Responses 
to these environmental variations are not available in the existing 
formulation of the circulatory system model. The implications of 
providing an indirect environmental stimulus to the circulatory 
system via the respiratory system enhances the capabilities of the 
integrated system. Preliminary simulations of this classification 
of experiemnts have indicated a need for additional modifications 
of the circulatory system so that it is sensitive to the environ- 
mental sti mul i . 

1 .4 Improved Respi ratory Frequency Expression 

Exercise presents the organism with the task of rapidly and 
optimally adapting to the requirements imposed upon it. One of 
these tasks is to increase ventilation in order to meet the ac- 
celerated metabolic demand for 0 2 by the exercising muscle systems 
and the corresponding need to vent the C0 2 byproduct. 

A ventilatory stimulus could be considered a physiological 
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variation which acts upon the respiratory centers and supplies 
information about the respiratory requirements of the body. There 
are two cl assi fi cati ons of ventilatory stimuli: humoral and 

neural. Humoral stimuli involve modifications in the physical or 
chemical properties of the circulating blood that stimulate the 
respiratory centers either directly through the blood or via 
afferent nerve endings located in receptors in contact with the 
vessels. These are termed chemoreceptors and baroreceptors . The 
other stimuli are classified as neural stimuli. They originate 


in the brain or in cutaneous, mucosal, or deep peripheral receptors. 
Generally, these respond to localized conditions. Humoral and 
neural factors -are postulated as controlling respiratory frequency. 

The complexities of the interactions between neural and humoral 
control are manifested during exercise. The individual must in- 
crease ventilation to a level sufficient to provide the additional 
0 2 intake and C0 2 venting demanded by the body's high level of 
exercise metabolism. However, exercising subjects ventilate at a 
rate that exceeds the increases that would-be dictated solely by 


arteri al 



pH, 


and P n alone. (7) 
u 2 


The difference between 


the humoral requi rements and the actual ventilatory rate must then 


be a function of neural stimuli. This could be the result- of neural 


pathways between the motor cortex and the respiratory centers and 
feedback from the foremen ti oned proprioceptors. The neural fac- 
tors assume even greater significance when ventilation is plotted 
as a function of exercise level and duration; however, there seems 
to be an upper limit of neural stimulation. It is noted-that there 
is a rapid increase in perspiration immediately after the onset of 
exercise; "this component has been postulated to be strictly neuro- 
logical in nature. 
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An important question about ventilation, both for resting 
states and during exercise, concerns the actual frequency of res- 
piration. It is obvious that identical ventilation requi remen ts 
can be fulfilled either with a series of rapid, shallow inspirations 
or with long, deep ventilations. 

The ventilation formulation in the respiratory system model 
is a satisfactory representation. The respiratory frequency ex- 
pression is the one for which modifications are proposed. It 
is desirable to'have a physi ol ogi cal ly based expression for res- 
piratory . freq uency since this could become an important inter- 
facing component for the respiratory-pulsati le cardiovascular system 
model . 
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2. EVALUATION OF EXERCISE SUBROUTINE MODIFICATIONS 
IN THE RESPIRATORY SYSTEM MODEL 

2 . 1 Statement of Objecti ves 

With the respiratory control system model functioning as an , 
independent system certain modifications of the exercise subroutine 
allow for a more realistic response to exercise stimulation. These 
modifications apply to both on- and off-transient responses. Eval- 
uation of these modifications which preserve the neural and humoral 
control of ventilation is desired. 

2.2 System Response Before Impl ementati on of Modi fi cati ons 
The discussion in this section refers to the respiratory 

system functioning independently with all modifications relating • 
to the computer program as listed in the Appendix of Gallagher's 
report. (7) Some of the modifications are included in the inte- 
grated ci rcu latory-respi ratory system model in addition to the 
modifications required for the interfacing of the two systems. 

For implementation and evaluation of exercise subroutine modi- 
fications for the respiratory control system simulations a- spectrum 
of exercise levels were simulated. Appropriate magnitudes and 
lengths of exercise levels were implemented. These ranged in mag- 
nitude from the resting state to a submaximal exercise level of 250 
watts- and in durations of 2.5 to 12.0 minutes depending upon the 
exercise increment. 

Although the simulation provides approximately 60 physio- 
logically oriented output variables the following variables were 

i 

closely monitored so as to determine the exercise simulation de- 
ficiencies . 

(a) inspired ventilation rate {Vj, 1/min) 



- 9 - 


(b) cerebrospinal fluid H + concentration ( c cSF(H + ) 5 
nanomoles/1 CSF) 

(c) arterial 0 2 tension (P a ( C0 mm H g ) 

(d) arterial C0 2 tension (P a ( 0 j » mm Hg) 

(e) tissue 0 2 metabolic rate (MR-j-^g y 1/min) 

(f) tissue C0 2 metabolic rate (MRj^gg y 1/min) 

(g) alveolar respiratory quotient (Alv RQ) 

(h) cardiac output (Q, 1/min) 

(i) brain blood flow (Q B , 1/min) 

Figures 2 and 3 illustrate the type of responses generated 
by the system- before the modifications were introduced. In general, 
the greater the exercise level the more pronounced are the vari- 
ables' transient responses. If one would superimpose a variable's 
response for a wide range of exercise levels a family of curves 
illustrating trends in the variable's response to exercise stim- 
ulation would be evident. After careful examination of the system's 
responses to each exercise level the 200 watt level was chosen as 
the base run fromwhich exercise subroutine modifications were 
evaluated. 


For the following discussion of the simulation deficiencies 
refer to Figure 3. Inspired ventilation, Vj , contains a reasonable 
neurological component as illustrated by the immediate on-set of 
ventilation when exercise is initiated. 


This is followed by a slowly rising Vj response, the humoral 
component. Sufficient ventilation is not available in order to 
properly regulate the arterial 0 2 and C0 2 tensions as demonstrated 
by these responses. The inadequate blow-off of C0 2 is reflected in 
the increased concentration of the H + ion in the CSF compartment. 
Since the CSF compartment's H + ion is an important regulator of 
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Figure 2. Selected variables from the respiratory system 

simulation of an exercise (work load) excitation 
of 150 watts prior to exercise subroutine modi- 
fi cati ons . 



nan omol es 
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Figure 3. Selected variables from the respiratory system 

simulation of an exercise (work load) excitation 
of 200 watts prior to exercise subroutine modi- 
fications. 
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ventilation one would expect to see an increase in the ventilation 
rate; however, this regulatory component is of minor importance in 
exercise simulations with normal environmental conditions, i.e. 
physiologically compatible atmospheric gaseous mixtures. 

Also evidenced from Figure 3 is the simulation's inadequate 
regulation of cerebral blood flow. This response correlates with 
an insufficient ventilation rate. 

The off-transi ent exercise response provides a ventilation 
rate that exceeds the necessary rate. In other words, the ven.til-’ 
ation rate causes an increase in arterial 0 2 tension and a decrease 
in arterial C0 2 tension which could be smoothed with proper modi- 
fication of the exercise subroutine.' In addition, the piece-wise 
linear off-transient Vj response is not physiologically justifiable. 

2.3 System Response After Implementati on of Modi f i ca ti ons 

Three basic modifications were implemented in the respiratory 
system program yielding improved simulations. These changes are 
reflected in the system responses for the respiratory control sys- 
tem independent of the circulatory system simulation. Consequently, 
these program changes are not contained ‘in the program listing in 
the Appendix. Variations of some of these expressions with the in- 
tegrated system are discussed in Section 3. 

In compliance with Astrand and Rodahl (9) a change in .the func- 
tional relationship between steady-state 0 2 requirement and exercise 
level (SS02W (WORK) ) was made. The functional discontinuities listed 
in the statements of subroutine SS02W(X) were removed yielding the 
rel ati onshi p 

SS02W( X) = ( X / 75 . ) + . 2 1 5* ( 75 . -X) /75 . , 0$X<75 

« -.072 + X/70 . , 75 250 

= 3.5, X>250 for X = WORK (2-1) 
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For a 200-watt exercise level this change increased the steady- 
state 0 2 requirement by 4.5% to 2.785 1/min and correspondingly 
increased steady-state ventilation approximately 6.4%. (10) Thus, 
this' subroutine dictates the steady-state level of required 0 2 for 
an exercise stimulus and can be easily altered to meet individual 
demands . 

Another modification which improved the ventilation response 
and consequently the regulatory aspects of the system involved 
the exponential functional expression that describes the on- and 
off-transient ventilation responses. VTIME is an expression in 
subroutine R C 1 2 which indirectly describes the dependency of ven- 
tilation upon magnitude and duration of exercise' levels. The 
VTIME expression which was applicable for both the on-transient 
and off-transient responses to exercise stimulation was originally 
programmed as 

VTIME = TCT* ( CXT-TI ME ON ) /9 . 2 . (2-2) 

Modification of this expression yielded for the on-transient 
case the expression 

VTIME = 1 . 1-1 .1*EXP(-TCT*(CXT-TIME0N)/1 .92) (2-3) 

and for the off-transient case 

VTIME = 1 .1-1 .l*EXP(-TCT*(CXT-TIME0N)/3.84) . (2-4) 

Here, • 

1/TCT ~ ti me -constant associated with the exponential 
functions related to exercise levels. 

CXT = simulated time, and 

TIMEON = time for initialization, of new exercise level. 



- 14 - 


Figure 4 illustrates the responses utilizing the foremen ti oned 
modifications. The modifications retain the neurological component 
of inspired ventilation, Vj, and allows steady-state ventilation to 
be approached much more swiftly with a slight positive derivative 
in the ventilation rate prevailing until exercise is terminated. 

The faster response in ventilation rate, which is more acceptable 
physiologically speaking, provides good regulation of the arterial 
Og and CO 2 tensions during the initial portion of the exercise 
stimulus in addition to good regulation of cerebral blood flow. 

Although not a critical problem, the off-transient response 

4. 

needs some refinement. The build-up of H concentration in the CSF 
compartment suggests further modifications in the metabolism form- 
ulations. Cerebral blood flow, although not unstable, possesses 
a response which should be smoother. 

To further illustrate the capabilities of the system to re- 
spond to a variety of exercise levels a 40-minute simulation in- 
volving the following series of exercise levels was run. 
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Figure 4. Selected variables from the respiratory system 
simulation of an exercise (work load) excitation 
utilizing exercise subroutine modifications. 




Time, min 0 • 5 10 15 20 25 30 35 

Work Load, joj 40 | 1 00 1 150 | 0 | 150 | 250 | 0 | 

watts 


Figure 5. Selected variables from the respiratory system simulation of several exercise 
(work Toad) excitations utilizing exercise subroutine modifications. 


<y\ 
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is regulated during on-transient periods and for low to medium 
exercise levels; however, the off-transient of Q B corresponding 
to sub-maximal exercise levels exhibits poor control which is sup- 
portive of the statements made in the previous paragraphs. 

The third modification which could be implemented without 
appreciable sacrifice in the simulation's fidelity involves the 
differential equation Subroutine RC13. Since the system's vari- 
ables are not rapidly varying as a function of time the thesis is 
the following. There is no justifiable requirement for a differ- 
ential equation subroutine having the capabilities of a 4th_ order 
Runge-Kutta method and an Adams-Moutl on predictor-corrector scheme. 
Using this criterion the following subroutine may be substituted 
for the foremen ti oned method. Note that the program listing in 
the Appendix does not contain this substitution. 

C ( 35) = C (35) + C ( 36 ) 

CALL RC14 
DO 1350 I = 1, M 
1350 C ( I ) = C(I) + C (36 ) *DC ( I ) 

RETURN 

END 

Here C ( 35 ) = time, 

C ( 36) = time increment, 

C(I) = system variable, and 

DC(I) = derivative of system variable. 

With preliminary evaluations of this modification there was 
demonstrated a small variation in the responses of the simulations 
using the two differential equation routines in the initial transient 
responses. Since the variables are changing most rapidly during this 
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time interval this phenomenon should be expected. With due con- 
sideration given to the variety of simulations that have been per- 
formed with the respiratory model it appears that this modification 
would present' the most difficulties when simulating extreme vari- 
ations in environmental gaseous concentrations. Simple exercise 
stimulus variations should be handled without difficulty. 
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3. INTEGRATION OF PHYSIOLOGICAL SYSTEMS 

3 . 1 Statemen t of Objecti yes 

■One of the major objectives of the overall research effort is 
to identify appropriate variables and subsystem components that can 
be utilized in an integration scheme for the physiological system 
models. Specific concentration within the present research effort 
is with the circulatory and respiratory systems. Supporting the 
objective is the establishment of an evaluation procedure. With 
exercise being the primary stimulus the proposed integration scheme 
for the ci rcul atory -respi ratory system is evaluated for both nor- 
mal and stressed physiological system states providing partial 
fulfillment of the objectives. Section 3.5 presents descriptive 
material aligned with the objectives of the integration plans in- 
volving other physiological systems. 

3.2 Ci rcul atory - Respi ratory System Inteqrati on Scheme 

In Section 1.3 a description, including a block diagram (Fig- 
ure 1], of the integrated ci rcul atory-respi ratory system model was 
presented. The three stimuli - .exercise, parameter, and environ- 
mental variations - were discussed with reference to their impor- 
tant forcing capabilities. This section contains a discussion of 
the interfacing component of the overall system. Refer to Figure 1 
for the pictorial representation. 

With regard to the exercise phenomenon it is important to con- 
sider both aerobic and anerobic oxygen deficits and debts. As an 
improvement' upon the formulation of the 0 2 metabolic rate for the 
tissue compartment- (humoral forcing component of inspired ventil- 
ation) in the respiratory system simulation the following interface 
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was established. The total 0 2 metabolic rate for the body is 
gi ven as 


(RM ro ~ oo DOB) = RMT < 2 ) + c ( 26 ) 


MR T(0 9 ) + MR B(0 9 ) 


( 3 — 1 ) 


whe re 


RMO 


DOB 


rate of 0 2 deli 


RMT{2 ) 
C ( 26 ) 


blood to the muscle 

tissue , ml /min . , 

rate of 0 2 delivery by the blood to the non- 
muscle tissue, ml/min. s 

M R t(0 j = metabolic rate of 0 2 in the tissue 
compartment, 1/min-. , and 

MR b( o ) = metabolic rate of 0 2 in the brain com- 
partment, 1/min. 


RMO and DOB taken from the circulatory system are functions 
of several physiological variables; thus yielding a more reasonable 
description of 0 2 requirements during exercise. The calculation 
of the metabolic production rate of C0 2 in the tissue compartment 
is retained in the respiratory system simulation. Also, direct 
neurological control of ventilation related to exercise levels 
is retained in the respiratory model. 

The total cardiac output is an important component of both 
systems. The circulatory system simulation describes cardiac out- 
put as a weighted expression of several physiological attributes. 

Here , 

QLO = (QLN) (LVM) (HSL) (AUH) (HMD} (HP L) (3-2) 

where 

QLO = left ventricle output, 1/min 

QLN = output of left ventricle under normal conditions, 

1 /mi n . , 
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LVM = effect of arterial pressure loading factor on 
left ventricle, 

HSL = basic strength of left ventricle, 

AUH = degree of autonomic stimulation of left ventricle, 
HMD = degree of deterioration of left ventricle caused 
by low coronary blood flow, and' 

HPL = degree of hypertrophy of left ventricle. 


Refer to Guyton (3) and White (6) for a detailed description of 
the terms in Eq . (3-2). 

In the respiratory system model cardiac output, Q = C(10), 
is described by a first-order differential equation and depends 
only on specific levels of arterial Og and CO,, tensions. (2) 

Based upon the comparison of these two formulations and the complex- 
ity of altering each, the decision was made to allow the respir- 
atory system component to receive cardiac output from the circula- 
tory system component. The formulation of cerebral blood flow is 
retained in the respiratory system. 

Another interfacing component involves the blood oxygen cap- 
acity. In the respiratory system this component is a constant 
{ H b = C(17)) alterable by an input data card. In the circulatory 
system this term is continuously calculated, thus it is provided 
to the respiratory system as 


where 


HM • OSA 
200 


Hb = C ( 1 7*) 


(3-3) 


HM = hematocrit 

OSA = arterial oxygen saturation, and 

Hb = arterial blood oxygen capacity, 1 0^/1 blood. 



In turn, the oxygen volume attached to hemoglobin in the aortic 
blood is developed in the respiratory system. The circulatory 
system utilizes this term in' the expression 
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OVA = 1000 • CHBA 


where 


CHBA - C a ( |^ b 0 ) “ ar ucri a i neiiiuy i vu i ii luiilsii u a u un , 

1 0 2 /l blood, 

OVA = oxygen volume attached to hemoglobin in aortic 
bl ood . 

This latter interface relationship is important in providing 
a possible pathway to the circulatory system for an environmental 
stimulus (alterations in gaseous • composi ti on ) . Since the venous 
blood in the respiratory model is not correlated with any partic- 
ular venous site in the ci rcul atory -model the calculations of 
tissue venous hemoglobin concentration )) is retained in 

the individual models. The foremepti oned features of the inter- 
face system contribute to a more physiologically reliable repre- 
sentation of the variables than formerly existed in the individual 


= arterial hemoglobin concentration. 


systems. 

The complete computer program for the integrated system is 
listed in the Appendix. However, there are some modifications 
which should be mentioned in comparing the individual respiratory 
system program with the respiratory system component of the inte- 
grated program. The original format and namelist statements are 
changed to comment statements. 

The following comment statements are listed at the beginning 
of the respiratory system program, SUBROUTINE GRODIN. These iden- 
tify the flow of variables between the two systems. 



- 23 - 


FOLLOWING FROM GUYTON TO GRODIN, 


GUYTON 

GRODIN 




QLO 

C( 10) 

{ URZ1 ) 

CARDIAC OUTPUT. 

( RMO + DOB ) /I 000 . 

RMT(2)+C(26) 

{ URZ2 ) 

TOTAL 

METABOLIC 



- 

RATE 

OF BODY. 

(HM*0SA)/200. 

C ( 1 7) 

( URZ3) 

BLOOD 

OXYGEN CAPACITY. 

URZ4 

URZ4 

(URZ4) 

FLG . 1 

= 1ST. TIME 


GRODIN CALLED. 

2 = GUYTON JUST 

OUTPUT. 

3 = GUYTON JUST 

OUTPUT AND 
DETECTED 
END RUN. 

0 = NOT ABOVE. 

URZ5 URZ5 (URZ5 ) WORK LEVEL (WATTS) 

FOLLOWING FROM GRODIN TO GUYTON 
GUYTON GRODIN 

OVA/lOO 0. CHBA (RUZl) OXYGEN VOLUME IN 

AORTIC BLOOD. 

In Subroutine RC12 the Guyton segment of the program is sent 
the work -load change. 

In order for the circulatory and respiratory systems' pro- 
grams to be compatible with their formulation of oxygen require- 
ments for various exercise levels, modifications were made in 
Subroutines FUNCTION SS02W(X) and FUNCTION SSVENT(X). The original 
Subroutine’ FUNCTION SS02W(X) is presented below. 
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FUNCTION SS02 W( X) 

C CALCULATION OF STEADY-STATE OXYGEN REQUIREMENTS FOR VARIOUS LEVELS 
C OF WORK LOAD (X = WATTS). 

I F ( X . GT.250. ) GO TO 1 
I F ( X . LT.75 .) GO TO 2 
SS02W = -.072 + X/70. 

RETURN 

1 SS02W = 3.5 
RETURN 

2 SS02W = (X/75.) + . 21 5* (75 . -X) /75 . 

RETURN 

END 

The modification of the above subroutine yields 
FUNCTION SS02W(X) 

C CALCUALTION OF STEADY-STATE OXYGEN REQUIREMENTS FOR VARIOUS LEVELS 
C OF WORK LOAD (X = WATTS) . 

SS02W = .195 + (X/84.15) 

IF (X.GT.210.) SS02W = 2.7 

RETURN 

END 

This modification produces a reduced level of steady-state 
oxygen requirement compared to the original respiratory program. 

The original Subroutine FUNCTION SSVENT(X) is presented below. 

FUNCTION SSVENT(X) 

C CALCULATION OF STEADY-STATE VENTILATION RATE AS A FUNCTION 
C OF TISSUE OXYGEN METABOLIC RATE 
IF (X.LE, .215) SS VENT « 5.398 



- 25 - 


IF((X.GT. .215) . AND.(X.LT.2.))SSVENT = 25. *X 
IF (X . GE . 2 . ) SS VENT = 50 . +50 . * { X-2 . ) 

RETURN 

END 

The modification of the above subroutine yields 
FUNCTION SSVENT (X) 

C CALCULATION OF STEADY-STATE VENTILATION RATE AS A FUNCTION 
C OF TISSUE OXYGEN METABOLIC RATE 
I F( X . LE . .195)SSVENT = 5.398 
IF (X . GE ,2 . ) SSVENT = 55.36 + 50.*(X-2.) 

IF((X.GT. .195) ,AND.(X.LT.2.))SSVENT = 27.68*X 

RETURN 

END 

A necessary condition that must be fulfilled for the integrated 
system to function properly is that both subsystems must be oper- 
ating with the same steady-state values prior to application of an 
exercise stimulus. After steady-state conditions are established 
for any particular type of normal or abnormal environmental or 
physiological condition, the integrated system components will 
function in concert for all parameter or stimuli disturbances. One 
key used in determining when steady-state conditions are reached 
is the observance of a minimal change in variable values C ( 1 ) - C(14) 
of the respiratory system model. These variables are related to 
the compartmental gaseous concentrations, cardiac output, and 
cerebral blood flow. 

Documentation of initial physiological' variable values are 
necessary so that both of the systems are functioning under com- 
patible initial conditions for all levels of exercise. As an 
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example of the type of documentation required, steady-state 
resting cardiac output is adjusted to 5.12 1/min in both models. 

All other variables which are dependent upon this variable are 
adjusted accordingly. Model modifications such as these are 
deemed necessary if both models are going to function in concert. 

The printed output data combines the forms of both individual 
systems. (6,7) In addition to the tabular form of the respiratory 
system model's output, the following circulatory system variables 
are printed. 

SECS = seconds 

PA = arterial pressure, mm Hg, 

QLO = left ventricular output, 1/min., 

PLA = left atrial pressure, mm Hg, 

PRA = right atrial pressure, mm Hg, 

VP = plasma vol ume , 1 , 

VPF = free fluid in interstitial spaces of lungs, 1, 

VTS = total interstitial fluid volume, 1, 

VUD = urinary output, 1/min., 

RMO = rate of oxygen delivery to muscle tissue, ml /min., and 
DOB = rate of oxygen delivery to non-muscle tissue, ml/min. 
Also, special plotting routines are available and provide for ex- 
cellent qualitative analysis and indicate variable trend setting 
phenomena. This feature of the simulation set-up is very useful 
when one is observing trends in a variable's response to a family of 
excitations, i.e. exercise levels. 

3 . 3 Exe rcj se Response for Nonstressed Phys i ol ogi cal System States 
As described in the previous section, in order to establish 
a base run, appropriate initial conditions have to be obtained for 
the respiratory system component. Basically, cardiac output and 
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metabolic rates must be compatible between the two systems. 

The input data as shown in Appendix 6.1, Table 2 of the refer- 
ence by Gallagher (7) was utilized with the following modifications. 



Input Data 




Card No. 

Va-ri able 

New Value 


10 

Q 

5.1200 


25 

MRB(C02) 

, .0450 


26 

MRB ( 02 ) 

.0450 


31 

FI (C02 ) 

.0004 


32 

FI ( 02) 

.2096 


33 

FI(N2) 

. 7900 


39 

PRINT AL TIM 

.2000 


45 

RMT ( C02 ) 

.1716 


46 

RMT( 02 ) 

.1950 


Table 1 . Input 

data cards reflecting changes in cardiac 


output and metabolic changes under 

normal environmental 


gaseous conditions. 



This simulation was 

allowed to run until 

steady -state conditions 

were 

reached, i.e. changes in C(l) - C ( 1 4 ) were minimal. These 

values for C(l) - C ( 1 4 ) were then established 

as the input data 

for 

normal environmental 

conditions. The variable values obtained 

for 

the steady-state conditions are listed in 

Table 2. 


Input Data 




Card No. 

Vari able 

Value 


1 

FA(C02) 

.1767 


2 

FA(02) 

.5338 


3 

FA(N2) 

.2895 


4 

CB (C02) 

.6345 


5 

CB ( 02 ) 

.0012 


6 

CB ( N2 ) 

.0011 


7 

CT(C02) 

.6142 
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Table 2 Continued 

Input Data 


Card No. 

Vari able 

Val ue 

8 

CT { 02 ) 

.0014 

9 

CT(N2) 

.0013 

10 

Q 

5.1554* 

1 1 

QB 

.7391 

12 

PCS F ( C02 ) 

46.3498 

13 

PCS F( 02 ) 

38.4441 

14 

PCSF (N2 ) 

70.6931 

Perhaps Q = 

5.12 should be used 

for complete 


bility with the circulatory system. 


Table 2. Initial condi t ions for physiological vari- 
ables under normal environmental conditions. 


Included here is a brief description of a simulation run 
for the nonstressed physiological system which correlated very 
well with a simulation- performed with the individual respiratory 
system. Refer to Figure 5. The stimuli were a series of exer- 
cise levels with specific durations as shown here. 


Ci rcul atory 
System 

Respi ratory 
System 


Exerci se 

Level 

Du ration 

(EXC) 

(Work load, watts) 

(min) 

1.000 

0 

0.5 

12.900 

40 

• 3.0 

38.500 

100 

3.0 

65.880 

150 

5.0 

1.000 

0 

9.0 


For each of the exercise levels and their corresponding trans- 
ient and steady-state responses the integrated system provided a 
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slightly increased cardiac output, inspired and expired ventilation 

rates, alveolar RQ, tissue P n and decreased arterial and tissue 

2 

arterial P n , and cerebral blood flow. All other physio- 


C0 2 u 2 
logical variables were altered accordingly. Although no definite 

trend was established, the VI-VE difference was altered when com- 
pared to the independently functioning respiratory system. 

Similar types of unstressed physiological system simulations 
were performed under altered environmental conditions. To eval- 
uate the initial conditions for the respiratory system model the 
input data as shown in Appendix 6.1, Table 2 of the reference by 
Gallagher (7) was used with the following modifications. 


Input Data 
Card No. 

Vari able 

New Value 

10 

Q 

5.1200 

25 

MRB(C02) 

.0450 

26 

MRB ( 02 ) 

.0450 

30 

B 

260.0000 

31 

FI ( C 02 ) 

.0192 

32 

FI ( 02 ) 

.7000 

33 

FI { N2 ) 

.2808 

39 

PRINT AL TIM 

.2000 

45 

RMT ( C02 ) 

.1716 

46 

RMT ( 02 ) 

.1950 


Table 3. Input data cards reflecting changes in cardiac 
output and metabolic changes under altered environmental 
condi ti ons . 

This simulation was allowed to run until steady-state con- 
ditions were established. Input data cards C ( 1 ) - C ( 1 4 ) shown 
below correspond to steady-state conditions for the altered en- 
vironmental conditions. 
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Input Data 
Card No. 

Variable 

Val ue 

1 

FA( C02 ) 

.1767 

2 

FA(02 ) 

.5338 

3 

FA{N2) 

.2895 

4 

CB ( C02 ) 

.6345 

5 

CB { 02 } 

.0012 

6 

CB (N2 ) 

.0011 

7 

CT ( C02 ) 

.6142 

8 

CT(02 ) 

.0014 

9 

CT(N2) 

.0014 

10 

Q 

5.1554 

1 1 

QB 

.7391 

12 

PC S F ( C 02 ) 

46.3498 

13 

PCS F( 02 ) 

38.4441 

14 

PCS F( N2 ) 

70.6931 


Table 4. Initial condi tions for physiological variables 
under altered environmental conditions. 

As illustrated by the values of Table 4 the combination of 
barometric pressure and volumetric gas fractions of the altered 
environment provided a physiological condition similar to that 
experienced at sea level. An interesting experiment involved 
implementing the conditions as illustrated in Tables 3 and 4 and 
observing the effects of exercise upon the system.. On- and off- 
transient exercise conditions were simulated for levels of 50, 

100, 150, and 200 watts. Each level was of a duration that allowed 
for the physiological variables to attain or approach their steady- 
state conditions. Results were similar to those obtained for the 
normal environmental conditions using the integrated system or 
for the respiratory system functioning independently. 
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3.4 Exerci se Response for Stressed Physi ol ogi cal System States 
A si gni f i can t - benefi t is achieved in using the integrated 
system as a means of determining respiratory system response to 
exercise when the circulatory system is not optimally functioning. 

In this way one can demonstrate the dependency of the respiratory 
system, in particular the ventilation rates, upon abnormal (stressed) 
functioning of the circulatory system. Thus, by observing easily 
monitored physiological variables noninvasive evaluation of cir- 
culatory system malfunctioning is achieved. 

The simulations involving exercise stimulation and the 
stressed circulatory system are implemented in the following manner. 
Under the influence of normal or physiologically compatible envir- 
onmental conditions the circulatory system is stressed (parameter 
variation). Examples of stressful situations include malfunctioning 
of the renal system, excessive fluid intake, unilateral o.r bilateral 
heart failure, regional volume loading, anemia, or hematocrit level 
variations. All of these are situations that could have been 
encountered in recent spaceflights or while performing normal 
routine earth-bound tasks; consequently, interest is keen on such 
simulations. Steady-state conditions are then established for the 
stressed condition with the appropriate variable changes realized 
for the respiratory system. The exercise stimulus is then applied 
to the integrated system. 

Justification for the above sequence of steps i.s stated in 
the following manner. Two types of systems are involved - a long- 
term model (circulatory) and a short-term model (respiratory). For 
a given stressed circulatory system state several hours or days of 
simulated time might be required for the system to approach steady- 
state conditions. Simulation with the integrated system is not 
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necessary for this period of time. Only after stressed steady- 
state conditions are reached is it necessary for the integrated 
system to become operative. Exercise stimulation is associated 
with rapid transient responses. It is applied for short periods 
of time; therefore, the short-term model is compatible for this 
time segment. 

Examples of variable responses for a stressed circulatory 
system under the influence of an exercise stimulus are shown in 
Figures 6-8. The stressed circulatory state was a bilateral 
heart failure. The pumping capabilities of both the left and 
right ventricles were reduced, simulating a mild heart attack. 

A simulated recovery period of approximately 2 weeks established' 
the stressed steady-state conditions. Responses to various levels 
of exercise were observed with a comparison of 100 and 150 watt 
exercise stimuli illustrated in Figures 6-8. -The exercise stim- 
ulus was applied for 5 minutes followed by a 5 minute period of 
no exercise. Included in these illustrations are only a few of 
the more significant variable responses. Direct comparison of 
the responses for the two exercise levels illustrates the depen- 
dency of each variable upon exercise stimulation for the particu- 
lar stressed state. A more detailed comparison may be observed 
from the tabular output data. Many details need to be refined 
before this particular system can be faithfully used in a quanti- 
tative predictive and/or diagnostic role. However, preliminary 
studies have been encouraging. 

3 . 5 Interfacing Possibi 1 ities for Physiological Systems 1 Models 
This section contains a brief summary of some possible inte- 
gration schemes for the respiratory, cardiovascular, thermoreg- 
ulatory, and the circulatory (long-term) system models. Even 
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though there remains considerable uncertainty involving the de- 
tails of the interfacing plans the overall goal hassome defini- 
tive aspects. The objective is a whole-body algorithm which sim- 
ulates physiological system responses to specific conditions and 
excitations directly related to those encountered in recent man- 
ned spaceflights (Skylab). 

Since there is a logical separation in the classification of 
models, long-term and short-term, it seems advisable to mold a 
total system which minimizes the interfacing structure and uses 
the aspect of simulated time to its best advantage. Thus, a 
plan which seems most encouraging at this time utilizes a separ- 
ation between short-term and long-term models. Each of the short- 
term models will interact with the other short-term models yielding 
appropriate transient responses for the stimuli. 

The respiratory and pulsatile cardiovascular system models' 
interface relies heavily upon compartmental blood flows and blood 
gas (O 2 and COg) concentrations. The respiratory model would re- 
ceive the following variables: cardiac output, cerebral blood 

flow, pulmonary vascular volume, and metabolic rate. The role of 
metabolic rate as a function of exercise would be handled in the 
respiratory system in a similar manner as described in the inter- 
face of the ci rcul atory-respi ratory system (Section 3.2). 

The respiratory model would be responsive to environmental 
conditions. It would supply the cardiovascular system with a 
physiologically justified description of respiratory frequency, 
arterial CO^ and O 2 concentrations, and an intrathoracic pressure 
variable. Some of the weaker features of the respiratory system 
model will be enhanced by inputs from the cardiovascular system. 

A similar statement applies for transferal of information in the 
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opposite direction. The two models must be made compatible for 
the resting steady-state case. To allow for this compatibility 
to exist it is necessary to establish the empirical a-v 0 2 dif- 
ference curve for a range of exercise levels as formulated by the 
individual respiratory program. 

A minimal degree of compartment redefining seems necessary 
with the above approach. If additional compartmental i zati on is 
deemed necessary the tissue compartment of the respiratory system 
could be subdivided to correspond to the legs, abdominal cavity, 
thoracic cavity, and upper extremities as described in the cardio- 
vascular system model. A restructuring of all blood flows and 
metabolic formulations would be necessary. With minor modifications 
the other two compartments, brain and lungs, and their associated 
variables could be related to the head and pulmonary segments in 
the cardiovascular system. 

The void of a compatible definition for compartmental blood 
flow and the dependency of compartments on the exercise phenomenon 
are two major obstacles in the integration of the respiratory and 
thermoregulatory systems. A more refined definition of blood 
flow may be achieved in the thermal model if a transition is made 
from lumped core regions to distributed core regions. However, 
this change would not yield itself to an easily implemented inter- 
face since all of the core, muscle, fat, and skin are lumped into 
a tissue compartment in the respiratory system. 

Under the influence of exercise it is reasonable to assume 
that the muscle compartment plays the dominating role as far as 
0 2 and C0 2 functions are concerned. Thus, the interrelationships 
among blood level parameters could be defined as follows. Arterial 
0 2 and C0 2 tensions could be supplied to the muscle model along 
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with a fractional supply of total cardiac output. In return the 
muscle compartment would supply a description of Og and COg 
metabolic rates, and 0 2 and COg gas tensions which would be inter- 
preted as tissue compartment variables in the respiratory system. 

The importance of the effects of the environmental temperature 
should be included in the respiratory component of the respiratory- 
thermoregul atory system. Temperature effects on total cardiac 
output (related to skin blood flow), ventilation rates, and meta- 
bolic rates are important. Additional terms or functions which 
are proportional to selected temperature deviations might be 
added to the calculation of cardiac output and the ventilatory 
controller equation to compensate for environmental temperature 
vari ati ons . 

Sections 3.2 - 3.4 presented an integration of the respiratory 
and circulatory systems. It allowed for a minimal reorganization 
of the existing systems and yet simulated with reasonable accuracy 
the influence of the exercise stimulus. Further compartmen tal i - 
zation of the respiratory system would no doubt be a part of the 
next level of complexity if the integration effort between these 
two systems were extended. With emphasis on the exercise stimulus 
the tissue compartment in the respiratory system would be sub- 
divided into (1) muscle, (2) renal, and (3) non-muscle, non-renal , 
and fat compartments. This partitioning would be in agreement with 
the major compartments of the circulatory system model. . Blood 
flow and metabolism for' each of the compartments would be developed 
as functions of Og requirements and COg production in addition to 
normal basal levels and neural control. The combined blood flow 
of these compartments plus the cerebral blood flow would total 
the cardiac output. In comparison to the interface system 
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described in Section 3, total cardiac output might be handled 
in the same manner. 
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4. RESPIRATORY FREQUENCY FORMULATION 

4 . 1 Statement of Objectives 

■Respiratory frequency is a term which can be related to the 
establishment of a sufficient Og supply and the proper venting 
of CO 2 . It is also an important driving force in the pulsatile 
cardiovascular system thus possessing the potential use as an 
interfacing component for the respiratory and pulsatile cardio- 
vascular system model. Consequently, a physi ol ogi cally compat- 
ible expression for respiratory frequency is desired. 

4.2 Prop osed Respi ratory F req uency Expressi on 

Separation of the factors causing ventilation rate, respir- 
atory frequency, and the different pulmonary lung volumes is 
practically impossible. Some humoral and neural control mechan- 
isms have been experimentally justified while many others remain 
undefined . 

Carbon dioxide can stimulate the respiratory centers dir- 
ectly and indirectly through the carotid and aortic chemoreceptors . 

For example, hyperventilation which decreases P/^ ^y several mil- 
limeters of Hg results in a decrease in ventilation. Variations in 
CO 2 concentration necessarily involve changes in blood pH due to 
the combination of CO 2 with HgO to form carbonic acid, and the 

4. 

subsequent dissociation to form H ions. This confounds efforts 
to distinguish CO 2 from the H + ion effects. There are small chemo- 
sensitive respiratory areas on the surface of each side. of the medulla 
oblongata near the entry point of the 8th, 9th and 10th cranial nerves 

4- 

This region is very sensitive to changes in H ion concentration 
in the cerebrospinal fluid (CSF). Because CO 2 diffuses from the 
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bloodstream into the CSF compartment much more rapidly than do 
free H + ions, the stimulus to the medullary centers is primarily 
due to C0 2 - The impact of the CSF compartment's H + ion concen- 
tration is discussed in detail by Gallagher (7) including changes 
in its weighting in the ventilation controller equation. 

Changes in arterial 0 2 concentration to levels below normal 
also have a strong stimulating effect on the chemorecep tors located 
bilaterally in the bifurcations of the common carotid arteries 
(carotid bodies) and along the arch of the aorta (aortic bodies). 
Their afferent neurons pass to the medulla through Hering's nerves 
to the glossopharyngeal nerves and through the vagi respectively. 

Understanding neural control of respiration is more elusive. 

In addition to neural pathways in the medulla that seem to have 
inherent inspiratory and expiratory patterns of rhythmicity, there 
are also other neural factors that influence respiration. Stretch 
receptors located in the visceral pleura and throughout the lungs 
give rise to the "Hering-Breuer " stretch reflexes whi,ch act through 
the vagus nerve to inhibit continued inspiration (or expiration) to 
certain limits of expansion or contraction. The Hering-Breuer 
reflex also contributes to cyclic respiratory rhythmicity. 

There are also other stretch receptors located throughout the 
body in joints and muscle tissues that give proprioceptive feedback. 
The immediate emotional state of the organism also acts to modify 
respiratory patterns. Speech, fear, rage and other stresses serve 
to alter respiration. 

One approach to the establishment of a respi ratory -frequency 
expression involves a work formulation. In this context work is 
associated with respiration itself. Otis, et al . (11) found that for 
low levels of ventilation the frequency of breathing could be fairly 
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accurately predicted by utilizing a minimum rate of work criterion. 
The problem definition assumed a sinusoidal airflow pattern, pas- 
sive expiration, a linear elastic element and nonlinear resistance. « 
Differential work was defined as 


dW = A ( t) + B ( t) + C ( t ) 

dW = KVdV + K ] a 2 sin 2 bt dt + K n a 3 sin 3 bt dt (4-1) 

where 

A ( t ) = elastic work 

B(t) = viscous work 

C(t) = work done in overcoming turbulent resistance. 


By graphing elastic, viscous, turbulent, and total work as 
functions of respiration frequency there emerged a frequency that 
minimized work for a constant alveolar ventilation. Also, it was 
determined that for too- low a frequency a greater amount of elastic 
work, was required to produce the large tidal volumes. In addition, 
for high frequencies, work was wasted in ventilating the dead 
space with each breath. 

Yamashiro and Grodins (12), assuming negligible turbulent 
resistance (C(t)-»0) and noting that the terms in the series expan- 
sion of B(t) were orthogonal with respect to the period T, minimized 
the rate of work done in breathing by using the equation 


dW _ 


Wf 


KfV T /2 + K 


Z 

1 = 1 


74 


(4-2) 


and the constant alveolar ventilation and dead space constraint, 


■V 


T 


+ V f - 


(4-3) 


Terminology used in Equations 4-2 and 4-3 follows: 
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f = respiratory frequency, bpm 

V j = tidal volume, 1 

1C 1 = total respiratory resistance 

■f* h 

a i = i 1 Fourier coefficient of airflow infinite series 
Vp = dead space volume, 1 

= minute alveolar ventilation, 1/min 

Utilizing an optimization procedure with Equations 4-2- and 
4-3 an optimal frequency was obtained as 

f opt = ((1 + 32 RCV a /V d ) 1/2 -1 )/16RC (4-4) 

with the same terminology as before in addition to 

RC = time constant associated with, resistance and 

elastance of the tissues. (Experimentally deter- 
mined to be .015 min) 

Yamashiro and Grodins (12) stated that since the airflow pattern was 
assumed periodic with no specific describing function the optimal 
work was achieved with a constant flow rate during inspiration and 
a decreasing exponential flow rate during expiration. This was 
associated with passive breathing. 

Examination of predicted optimal respiratory frequencies for 
optimal and sinusoidal airflow patterns indicate that sinusoidal 
airflow can require up to a 11.7% increase in work cost when com- 
pared to the optimal airflows. The situation is altered if active 
expiration is considered.. For active expiration the optimal air- 
flow pattern assumes the shape of a square wave with maximum and 

O 

minimum values of +2(V p f + V^) . Comparison with sinusoidal airflow 
indicates a 23.4% work cost increase if the pattern is sinusoidal. 

•Research involving active expiration seems to provide the 
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transition to the study of airflow patterns during exercise. Yama- 
shiro and Grodins (13) examined the minimization of work during ex- 
ercise when functional residual capacity (FRC), the airflow pat- 
tern, and respiratory frequency were included as controlled vari- 
ables. The result of their derivation yielded the optimal fre- 
quency expression 

f opt = (0+32(<l + a)/a)RCV A /V D ) 1/2 -l)/16((l + a)/a)RC. (4-5) 

The only term which has not been previously defined is "a". Here 
"a" is a ratio of inspiratory elastance (Kj) to expiratory elastance 

(k e ). 

When correlated with experimentally determined curves, a time 
constant (RC) of 0.015 min. was obtained. This compared with a 
similar value of 0.011 min. in research reported by Mead. (14) 

The value of "a" was assumed to be 1.95 based on data taken by Rahn, 
et al . (15) 

At levels of ventilation near rest, the experimental data 
appears consi stent 'wi th the expression for respiratory frequency 
associated with a sinusoidal airflow pattern and constant expiratory 
reserve volume (ERV). This is Equation 4-4. However, at greater 
respiratory rates, the expression which accounts for FRC (Equation 
4-5) provides the best fit with experimental data. 

Based upon the previous discussion and the importance of exer- 
cise in the simulations associated with the project's goals a 
recommendation to implement Equation 4-5 in the respiratory pro- 
gram appears reasonable. This formulation would be substituted in 
Subroutine RC12 for the expression 


FREQ = 8.1 + 7.815‘*(RMT(2) + C(26)) 


(4-6) 
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which gives an empirical representation of respiratory frequency 
as a function of 0 2 utilization. This change would require a modi- 
fication in the calculation of dead space ven ti 1 a.ti on . Using the 
discussions of the cited references an empirical formulation for 
dead space ventilation as a function of exercise level could be 
established. Upon implementation of the -appropri ate variable desig 
nations an improvement in the calculation of respiratory frequency 
is proposed. 
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5. CONCLUSIONS AND RECOMMENDATIONS 

Modifications of the individual respiratory system model have 
produced a system which simulates the effects of exercise in a 
manner which is acceptable with regard to regulation of Og and CO 2 . 
The neural and humoral control of ventilation during on- and off- 
transient exercise stimulation is justified. Likewise, other vari- 
ables are altered according to the regulation of gas tensions and 
ventilation rates. 

Improvement of physiological system models is always possible. 
Such is the case with the respiratory frequency expression. With 
projected plans for using the respiratory frequency as an inter- 
facing variable between the cardi 0 vascul ar and respiratory systems 
it is desirable that the expression be valid for the types of 
simulations performed. The proposed expression has met with success 
with various levels of exercise as discussed in the references of 
Section 4. However, there is no sound documentati on as to the per- 
centage contribution of humoral and neural control as given by this 
frequency expression. This constitutes another level of complexity 
in both experimental and theoretical respiratory research. 

To fulfill the goals of the research effort, future work should 
emphasize the integration of the short-term and long-term models 
so that an adequate whole-body algorithm is realized. There are 
many problems that will arise as the integrations of the inter- 
acting physiological systems are achieved. 

Results of this study are encouraging since a minimal amount 
of individual system modification was necessary for the integration 
of the circulatory and respiratory systems. Perhaps features of 
this approach can be applied to other integration endeavors. 
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Continued improvement of the integrated systems will help establish 
a better understanding of the interrelated physiological processes, 
especially those which are influenced by altered environments or 
stressed physiological states. Obviously, the continued develop- 
ment of appropriate data integration and display routines should 
parallel the modeling phases. 

The next phase of the project will undoubtedly be an enlighten 
ing experience with significant gains made in the establishment of 
techniques and useful procedures for integrating complex physio- 
logical systems. As an additional benefit of the modeling effort 
it is hoped that correlation of simulated and experimental data 
will provide potential noninvasive diagnostic and therapeutic cap- 
abilities for both stressed and nonstressed physiological states. 



6. APPENDIX 



C PROGRAM GUYTON 

G C t PCULAT fJRY DYNAMICS - CIPCLE 

C CIRCFl 

INTEGER URZ4 

PEAL LVM.I, IFP.LPO.KF.KEl.KOU.KIR.KU K I , KCO, KEO , KN t, XN 3 ' 

RCAL NAE.NEO.NIO.NOU.Il , LPK ,K II), MC2 ,M.Z ,KC Z, HPL ,H°R , 1 2, 1 l.R^n 
DIMENSION FUN LI14 ), FUN 21 LA) .FUNK 14) ,1 UNA ( 1A ) .FUNGI l4),rUN7tl4 J 
COMMON /ARR A Y/T . I ■ VBO , WS ,VP A, VAS , VL A . t RA , V AF , P A , °AM , L VM, 

* VRF.PRA.QRN.VPE .PPA.PPI , I.PA, RPA , RVM, VL E. PI A.QLN.Pl l , 

* ' Vld.RPV.RPT ,POL,«PO,«5 . VVE , VVB ,P VS. P3V ,P VG , OVC , AVF 

COMMON/ARRAY/CN2.CN3.RVS ,PGS , RT P, CAC.HRQ, Ql H, CVS , OP A, OAS ,OL A .OR 4, 

* PA1,AUC.4UP,AUN,AUA,AU2.i.U 3, OAU.AUJ.AU , AtJQ , AUH ,VVA , 

* AUR, AUM, AUA.VIF.PUl.PT T,1 TS.PIF.CP I .PTC, C“P,PPC .PVG 
COMMON /ARRAY/PC .PCO.VTC ,PLD , VTL . VT 0 PO, DPL , CP 1 , 1JPC, OP 1 ,1-PO.nLP , 

* OPP,CHY,P°M,PG.«. CPG.PGP.t F l ,PGX , PGC , PGH.PGZ , VGU , VG , 

» FPF,GP1,GP2,GP0,AA°,RR .' FN, APO.GLP.PFL.r.FP.TPR.VUD 

COMMON /ARR A Y/RFK .NOn.NEO.NAE.VFC.CKE.I 110, KFt , K l R,K 1 r, KCO .KEO.CKI , 

* CNS.CCD.VIO.SE ,NI ,VIC, l .VTY.Z , VT7 , VUZ , TVZ , FPZ , 

* OFZ.X ,12 ,PR t , VTS , VP . MRP, I FP , GPR, KN3.KN 1 , AMR , AMP 
C0MM0N/ARRAY/4 M l, AMC.AM2,AM3,AH5,AM .t.NE , AGK , ANP, AM , ANC , AN2 , A'13, 

* AN5.ANM.VB .HMl.MM ,1U I E , V I B, V IM ,RC 2, P02 ,RKC ,R C l , 

* P.CO,VPC,RSN,OVA,eFN,RC0, iUM, PIU , QS V . POT , POO , POP , AR I 

COMMON/ ARR AY/ AR2 . PCC , AP 3, ARM, CNH , GFN.-.H7 , A HR .AH .AHC.AH1.AH2 ,AHA,' 

9 AHM,CNY,CNX,VV1,VV2 ,VV5,"V6,VV7 .TVO.VTH.hSR.HSL.NIO, 

* -SR .VVR.RAR.CV .GN7.AUX, lUK.AUZ.Y , CFC , CPK , PLE ,CPR 
COMMON/ ARR AY/ LPK.OPLl.HYL ,KT0, AM T , ANT .t'UK, PON , A IK , 4 2K, A3X ,CNH ,CNZ, 

* AHK.SRK.VS .V20.Z1 , Z2 Ml ,ZA ,Z5 ,Z6 , Z7 ,Z9 ,HMK, 

* HXM.PIJV.PUZ ,R00 .OO2.R0F, «02 , POA , POY, AN1J , POR , GF2 .HMD 

CQMMON/ARRAY/OHP.POO, 13 ,U .VPl.Tl « »F3, GFA , AUP, A UV, R V 1 , AUY ,niJT , 

* OSP .AHZ.AHY.OSA.PPl ,CPN, 'US.PLF.PPO.PPN.POO.PFl .orp, 

* VPF, PPR.PMC.PMS, PMP.HR , IPF.PCP.OAl.OLZ.OPY.OPZ.GPZ 

COMMON/ARRAY/NOZ .KCZ.VIZ ,HPR .HPL.STH 1L 0, EXC ,02M, p A2, PP2, S VU , AUL , 

* VV9.O2A.0l .EXE ,ARF .ORF. (SM,i(FM,RAm,|)VS,PVC,1'M0,0CM, 

* PHC:,P2G,HF9,P00,P0E,AMM \AK,P0M,nMM,PM!,PH3,PMA,FXl 

C0MM0N/ARRAY/Q2 ,m ,PP5,PK1,Z9 ,Z10. II l,Z12, Z13.Z14, Z15.Z16.PK2, 

* PK3.FIS,STA,PAR,GRL.ANY, VNZ, ANX, ANV.ANN.ANR.AUO.AURt 

* AUS.A37B.H1 ,A2 , A3 .OU'HY U9 > , T ITL E I 400 ) , OlIMN Y( 40 ) 
COMMON/NUMERO/K.NO 120) ,NTIMEC , UNI TS ,M ,NT 1 MfcP , NN, MAXNC, NT [Mf 

1 , 1 PNEXT 

COMMON/STCRE/NGt .NG2.NG3 ,NG4 , NG5 . NG6 . \G 7 , NGR ,NG9 , OT , TLP .TNP.ND, 

* TM.TMM.NPIHST.ZZI 15) ,0L 'IV) ,(lUY<9) .YHIN(IO) .YMAXtlO), 

* N.PTI U),DETAaO),NGHAP-mO) .GRAPH) 10) .HrACI 19) .NOEXP 

* , 0 TMA X 
CCMMi^N/TAPE/TOTAL 

OATA FIJNl(l) .FU'll (2) , FUN It 3 ) , FIPJ ( ( A ) . -UNI ( 51 , TUN 1 ( 6 ) , FIJN1 1 7) , 
*F'JNlt « I, FUN 1(9 I , C UN1( 1C) .FUN 111 1 1 ,FU 111 12) ,FUM I 1 3 ) .FIINl ( 1 A ) / 

»0. ,1.0 A, AO., 1.325, 1 25... 9 7, 160. , . '13 , / 79. , . 59, 240. , 0. , 2 A 9. , 0. / 

DATA njM2< 1) , F UN2 ( 2 ) .FUN2 13 > » FUN? (A ) , : UN? 1 5) , FIIM2 1 ft ) , FUNZ (7 ) , 
*FUN2(8),FUN2( ) ),FI)N2( 10 ) , FUN 2 11 1 ) ,F u.t/ ( 12) ,FUN2( 1 3 ) ,F IJN? ( 1 A ) / 
♦—ICO. ,0.0, -fa. , 0. 0 ,-3. ..75,— 1. ,? .6, 2*. l.B. H., t 3.5, 1000.. 1 1.5/ , 

OATA FUN3 ( l ) , FUN 3 ( 2 1 , FUN 3 ( 3 1 , FIJN3 ( A ) ,PUN3(A) .FUN 3 ( ft ) , FUN3 ( 7) , 

* FUN3I R),FUN3t9),FlJN3(10),FUN3l 11), FUJI 121.FUN3I 13) ,FUN3t 1A)/ 
*0.0,1.06,20. ,.97,24. ,.93 ,30. ,.R ,3P. , . Aft ,45., 0. , A 5., 0./ 

OATA FIJNA ( 1 ),F1JN4(2J , FUNA I 3 ) , FUNAI A) ,=UN4(SI , FUNA ( 6) , PUNA ( 7) , 

* FUNAI i3 > ,FUNA [ 9 ) ,FUNA 110 ) , FUNA (11), FU'iAl^l.FUNAIU), FUN 4(141/ 

*-100. ,0. ,— A. ,0# ,-l. ,3. 6, 3,. 9. A, 6., 1 1 , ft , 10 . ,1 3 .5 , 1000. . 1 3 .5/ 

DATA FUN6 ( 1 ) ,FUN6( 2 ) . FUN61 3 3 , FIIN6 1 4 ) • FUN 6 1 5) .FUN6I 6) ,FUN6( 7 ) 

JFUNftf 3) .FUN619) »FUN6(10) ,FUN6 (11) ,FUNfr( 12),FUN6( 13 ) ,FUN6( 14 1 / 




-P" 

SO 

I 



*-100., IOO'JO.,0.,70., .4,9.3, .8| 3.3, 1. 1. 3, 1. 6, . 43 , 100. ,0./ 

OAT A FUN7 ( 1 ) , FIJN7 1 2 ) , FUN7 ( 3 ) , FUN7 l A ) , = UN7 ( 5) , FUN7 1 6) , FUN7I 7) , 
♦FUN7I Bl,FUN7(9l,FUN7( 10) ,FUN7(1 U . FUN U12 ) ,FUN7 ( 1 3 ) , FUN 7 (14)/ 
*0# 1 7 • f 30 * 1 6 » 2 5 » 60* » 3 • , 1 00 « 1 !•* 1 60 . , « l)fA00tf*0bt 400 • , • 05/ 


<50 CALL PUTIN 
U n Z4 = l 
C 

IFU .or. 0.5) 1*0.5 

c 


100 CALL PUTOUTIU«Z4,URZ5) 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FOLLOWING FROM GUYTON TO GROOIN... 


GUYTCN 

OLO 

[R*0*D0B 1/1000. 
<HM*OSA)/200. 
U»Z4 


UR Z 5 


GRCDtN 

CtlO) (UfU.)CARQIAC OUTPUT. 

RPT (2)+C(26) (URZ2IT0TAL METABOLIC RATE QF BODY. 


C ( 17) 
URZ4 


URZ5 


IUPZ3 ) 3L000 OXYGFN CAPACITY. 
(URZ4'=LG 1 = IST.TIHE GR001N CALLED 
2= GUYTON JUST OUTPUT. 
3=GUYT0N JUST OUTPUT 
ANT) DETECTED END RUN. 
0=NOT ABOVE • 

CURZ5M0RK LEVEL ( WATTS ) 


FOLLOWING FRCM GRCCIN TO GUYTON 

GUYTON GRODIN 

OVA CF8A4100Q. (RUZ1) OXYGEN VOLUME IN AORTIC 8L000 


C 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


102 UP.Zl = OLC 

(JRZ2 = IRMO ♦ DOR ) / 1000. 

URZ3 = (HM » OS A )/200 • 

CALL GROO INtURZl. URZ2.URZ3 tURZ4*URZ5. TUZ1 ) 
URZ4 3 0 

T = T ♦ 12 


CALL 

CALL 

4 

4 

• 

4 

4 

4 

4 


HE MG (FUNL.FUNJ ,FUN3 triJNA I 

FEMO IAM1,ANV,,ANU.ANY,ANZ , ARM , AUH , AU* , AUY, AVE , HF M ,BFN , 

CN2.CN3.CN7.CV f nAS,OLA t CPA,DRA,OVS,F IS.HMO.HPLi 
HPR ,HSL,HSR,I2 ,LVK,PA , PAM, PAP, PC , PGL , PGS ,P LA, 
PPA,»Pi,PP2,PPA,PP 1,4VS,OAO,OLN,OLO,QPC,ORF ,ORN, 
ORC,OVO,RAM,RAR,RBF.RPA,RPT ,RPV,RSM, RSN.RVG.RVM, 
RVS.U ,VAE,VAS,VBD , V IM, VIA , VLE , VP , VP t , VPE ,VK A, 
VRC,VPE,VVE,VVB,VVS,VV7,VV8,X , F UNI , FUN2 , FUN3 , 
FUN4 ) 


CALL AUTO 
4 
4 
4 


I AU ,AUR*AUC,AUM.AUJ, AUK, AULi AUM, AUN» AUD, AUP ,AUO, 
AUR,AUS,AUV,AUX , AUZ. AU4, AU&, AU8 , AL B, OAUtEXC ,EXE, 
FX l , 1 2 , PA ,PAl f PQ0*»OT,P20»STA,VVR,VV9,Y , 1 , 

ZB ,ZI2) 


If I I3.LE.I2IG0 TO 168 
IF ( AB S ( OAU-AUJ > .GT.CAl I GC TO 100 
UO IF (ABS(QAU-OLO).GT..2)GO TO 100 
IF UBS(CAC-QPO).GT,.2)GO TO 100 


1 

>01 

O 



IF (A9S(QA[)-0R0>.GT..4)GC Til 100 
C 

169 CALL HCRNCN . (AH , 4»C , AMP , AMR , AMT , All , AN", CKE.PA, Z 1 FUN7, 

* AGK,ANC,ANP,ANR,AN1 AN V , AMW , AN1 ,CNA, CNE ,CFN , 

* I ,RFK,A2 ,A3! 

C 

CALL BLOOD IHK’I.HM ,HMK,I ,POT POY , POl , P02 ,RCt ,RC2 , RCO, RKC, 

« V6 ,VIB,Vlt,VIM,VP .VRC1 

c 

CALL MUSCLE (ALOiAMMfADHt AU D , AAK BFM r EXC ,HM ,1 , MHO, PM"! ,OS A, 

* nVA,QVS,(l2A,PP0,PKl. PK 2i PK 3 , PMO, PM t , PH 3, PMA ,PM5, 

* 1 POE,PC1M,PVO,P?0,DOK RMO t VPF,Z5 ,Z6,RUZ1) 

c 

C4LL AUTQRG ( ACM , ARM , AR 1 , AR 2 , AR • . A IK, A2K , A3K.BFN , DOB , HM it, 

* MO 2 ,OSY ,OVA , Q2H,P07 P08, PQC , PODiPOK, PON, POR ,POT , 

* POV , POZ , PlOi Q02 tRDC, Z ,ZA , Z 7 > 

c 

CALL AOH (AH ,AMC,AHK,AHM,AHY AHZ , AH7 . AH8 , AU», CNA, CN8 ,CNR, 

* CNZ , I ,PPA,Z) 

C 

CALL MISC1 IAHM,AU4,AU6,I , SR SRK,STH,TVO,TVZ, VEC,VIC,VTW, 

* vvf, vv6 , W7, z ,vg> 

C 

CALL HEART ( AUR ,DHM , HMD, HR ,1 PA , PMC , PMP,PMS , POT, PRAiOAO , 

* OLO,RTP,SVO, VAE,VL( .VPE, VRE, WE, HU 

c 

CALL CAPMBD (BFN,CFC,CPI,CPP,DFP,I ,IFP,PC ,PCD,P! F,PL0,PPC, 
=> PRP,PTC,PTS,PTt,PVO.PVS,RVS,Tvn,VG ,VIO,VIF,VP, 

* VPD ,VTC ,VTD, VTL,VTS,VUD ,Z ,Z1 ,FUN6l 
C 

( = 1*1 .2+T-T1 

1I=ABS(VP1ZVP0/1I 

IFdl.LT.il I=!l 

ir t 13+T-T1.LT.1) I =1 3+T-T 1 

T=1*T1 

Tl=T 

c 

CALL PULMOM (CPF ,CPP ,CPN , DFP , I ,PCP,PF! ,Pl A ,PLF , POS , PP A,PPC, 

* PPfl.PPl ,PPN, PPO.PPR.VP , VPD i VPF ,Z ,Z31 
C 

CALL MI SC 2 (HPL,HPR,HSL,HSR,I,l*(,PPA,POT,STH,ZlO,ZU,Z13» 

C 

CALL PROTEN (CHY,CPG,CP( ,CPX,CPP,CPR,CPt,OLP,OLZ,OPC,nPI,DPL, 

* DP0,0PY,GPD,GPR,I ,IFP,LPK,PC ,PCE ,PGX , PRP , VG , 

* VTL ,Z i PPD1 

t 

CALL KIDNEY ( AAR , AHM , AM , API), ARC , AUM.CNE ,CN X ,C NY,GBL ,GF A ,GFR , 

* GF2 ,r,F3,GFA,GLP, t iNAF., NFI),N I D, NOO, NO Z ,P A ,PAR, 

* PFL,PPC,RBF,«EK,KFN,Rft , STH ,TRR , VI M, VUC , Z I 
C 

CALL ICKS (AM ,CCC,CKE,CKI,CNA,I iKCO.KE ,KE0,KI ,KI0,KIE, 

* K!R,Kni),NAE,REK,VEi:,VIC,VID,VP ,VPF,VTS,ZI 

c 

CALL GELFLO ( CHY ,CPG ,CP( ,GPR ,HYL, I FP, PGC , PGH.PGP, PGR ,PGX,PI F, 

* PRM.PTC.PTS.PTT.VG , VGD, V IF , VRS.VTS, V2 0, FUN6 1 
C 

GO TO ICO 
END 


I 

\J1 

H 

i 


SUBROUTINE PUTIN 



c 


COMMON/ ARRAY/AISOO) , TITLE IAOO) ,CDL(7l') , ALPHA! 20} 
CUM<1N/NU.MFR0/K,N0(20) .NTIMEC ,UNI TS, HZ , NT I MEP ,NN,MAXNO, NTIME 
I .1PNPXT 

C0“H0N/STCRE/NGl,NG2,NG3rNGA,NG5 t NG<. NG7 ,NG8 ,NG9 1 DT.TLP.TNP ,N0, 

* TM,TMM,NFIRST,ZZ< 15) ,0L Vt 9) , OBY t9) , YMINtl 0> .YMAXUO) , 

* N.PTtlS) ,BETAt 10) rNGRAPHUO) jGRAPHt 10) , HEAD! 19) ,NOEXP 

* .dtmax 

CATA All/'ALL * / , BLANK/ • •/, SAME/' SAME • / 

DO L J=l% AOO 

1 TITLEI JMHLANK 
KZ=0 

NT JMEP=1 
KN = 1 
PAXNC=l 

2 RFA0(5,10C) VAL UE, NUM3R0 , SYMBOL 

100 FCIR“AT <F13.6,2X, I5»2X,AA> 

IT ( 3AXN). LT.Nll Mt 5R0) HAXNC=NUMBRO 
1 F ( 'IJM«JRO . FO .0 ) GO TO 

At MW “SO I =VALUF 
T I TL S t N'JMA’O) =SY M ROL 
GO TO 2 

3 D5P=A ( 286) 

SEAntianil NOFXP, IHEA0(J|,J=1,19) 

101 FU»MAT(1A,19AA) 

36 ReAO(5t200) (ALPHAtJJ , J=l,20> 

200 FORMAT 1 20AA > 

201 IF ( ALPHAU ) -EO.SAME) GO TO 32 
IFlALPHAtD.NE.ALL) GO TO A 
RE AC ( 5,300 1 NT 1MFC, UNITS 

WRIT* (6,102) .N0EXP,(HEA0tJ),J=l,19) 

HR1TF(6,71I UNITS, ( TI TLEC4) , At J) , J = l , MAXNO) 

GO TO 31 
A 00 5 K*l,20 

IFt ALPHA (K).EQ. BLANK) GO TO 6 

5 CONTINUE 
K=2 1 

6 K=K-1 

00 10 J=1 ,K 
L= l 

7 1 F I ALPHA t Jl.EQ.TITLEU)) GO TO 9 
L=LH 

IF(L.LT.MAXNOtl) GO TO 7 
9 COL(J)=A(L> 

K0( J) =L 
10 COMT 1NIJE 
Gtl TO 3A 

32 00 33 4=1,K 
aleha(j) = reta(J1 
tiOt J1 sACiMAPHt J ) 

33 cru J)=r,PA5HiJi 

3A RE AO [ 5, 30C) NT IHCC, UNI TS,IPNEXT 
300 FORMAT) [6, A A, 12 1 

WRITE (6, 10?) NOEXP, (HEAOt J) ,J=1,19) 

102 F0RMAT(UU///2OX,«EXP * , IA, 19AA///1 
IF(K.GT.IO) GO TO 70 
WRITC(6,211 UNITS, (ALPHA!.)) ,J=1,K) 

21 FORi'iTt'O ',AA,10t6X,AA,lX>) 

WP1TFU.22) tCOL(J) ,J«l,K! 


I 

VJ1 

ro 

i 



22 FTRMAT ( • «,5X,'0'»2X,F10.4,9( lX.FlO.d) 

G: TO 31 

70 WTITE(6,71I UNITS, (ALPHA(J),CaLUI,,W,K) 

71 FORMAT ( 60X 1 2 HO , A4//5 1 4X, 44, • = • ,F 10 .4,4X> 1 
31 RETURN 

END 

SUBROUTINE PUTrUT(UKZ 4 ,URZ 5 ) 

INTEGER UPZ4 
C 

COHMCN/APRAY/A(400) , T1 TLE (400 ) ,COL (20) , ALPHAC20 ) 
CUMMON/NUMERO/K.N'MZCI , 

* NTI«EC.UNITS,NZ,NTIHEF*NN,HAXHC,NT1HE 
l tIPNEXT 

COMHON/STORG/NGl,NG?,«G3,NG4,NGS,NGZ,NG7,NGa,NG9,OT,TLP,TNP,NO, 

* TM,TMM,KF 1RST , ZZ ( 1 5 ) t QLY (9),0BY(9),YMIN( 10) , YM AXt 10) , 

* N,PT( 18) .BETAdO) .NGPAP'r (10) , GRAPH (10) ,HEACU9> ,NOEXP 

* i QTMAX 
COVMDN/TAPE/TOTAL 

OAT A S ECS/' SFCS*/,TM IN/ 'MINS' /.HOUR / 'HGUft* / .DAYS/' DAY S' / 

DATA ALL/ 1 ALL '/, PLANK/* •/ 

C 

C WATEXC(13,2) IS RELATION BETWEEN WATTS AND PARAMETER EXC (FRCM 

C OR. WHITE) TO ALLOW THE COMBINING OF ( RODIN WITH GUYTON. 

C 

DIMENSION WATEXCI13.2) 

DATA HATEXC/1.,5. , 10.,ZO.,3Q.,40.,50. ,60. ,70. , 80. , 90. , 10O., 1 20. , 

1 0., 16. 6, 32. 2,59. 1,93. ,103., 123., 14C. , 157. , 170. , 182. , 193. ,212./ 

C 

T*All) 

C 

C ' ICCNVT IS FLG l = CQNVERT EXC ( A f 3 9)) TO WATTS FOR GROOIN. 

C 0=00 NOT CONVERT EIIC TO WATTS. 

1C0NVT = 0 

IFMJRZ4 .Ft). I ) ICPNVT * l 

c 

NT I Pf = T/ 1440. 

IFtUNJTS .CO. SFCS) NT IMF = T 6 60. 

IFtUNITS .ft). TMIN) NTIME a T 
IFtUNITS .FC. HOUR) NTIME » T/60. 

IFtNTIME .IT. NTIMFP) GO TO 65 
C 

C HEBF IF IC PRINT. 

6 IF (UR Z 4 .NE . 1) URZ4 * 2 

I F ( ALPHA ( l | .ME. ALL) GO TO 7 

WR I TF (6,71 ) NTIME, (IN ITS, (TITLE(J),A( ) ) , J-l, MAXNO) 

GO TO 51 

7 CO 20 I « t ,K 

II * NC(l) 

COLtll . A( ui 

■ 20 CONTINUE 

(F(K .GT. 10) GO TO 70 

WP IT E ( 6 , 21 ) UNITS, ( ALPHA(J) , 0 = 1 , K) 

21 FOR V AT( ' 0 * ,AA,10(6X,A4,IX)I 

WRITE( 6 , 31 ) NTIME, (CCUJ) ,J= 1 ,K) 

31 FORMAT ( * M6,2X,FL0.4,9(1X,F10.4)> 

GO Tr) si 

70 WRITE(6,71) NTIMF, UNITS, (ALPHAt J) ,CCL1 J) , J*1,K) 

71 FORMAT ( / / /56X> 15, IX, A4// 5(4X , A4 , 1 * *,F10.4,4X)> 

Si NTIMEP => NTIME + IPNEXT 


t 

vn 

U> 


l 



c see if tine tc stop prfsfnt time stfp. 

53 I F I NT 1 1C..LT .NT IMEC ) GO TO 65 

5A REA1M5.A00) NTIMEC.CUNITS. IPNEXT, SYM3U.CVALUE 
AOD F0RMAT(16,AA,I2,AA,E13.6I 
C 

C BLANK T!«E STEP CAPD FNDS RUN. 

IF1 SYMROL.EQ.CUNITSI GO TO 66 
tFICUMTS.NE. BLANK) GO TO 59 
1FIM2J.GT..5) A(2>=.5 
A 50 no 55 MN=l,‘IAXNO 

I F ( SY M3GL .FC.TITLEIRM) GO TO 57 
55 CONTINUE 

57 W“ITE16,5R) N T IKE »UN I TS* SYMBOL »M MN) »C VALUE 

53 FI’RPATl/lSXt'AT'tlSflXtAAt' INTO THE SIMULATION, THE VALUE OF •, 
* AA.'WAS CHANGED FROM ' iFIO.3,' TO '^10,3/1 

A ( MM =CVALUE 

C SET FLO. INDICAT ION CHC-ED.WORK LOAD IF EXC INPUT. 

IFt-M .EO. 319) I CON VT = l 
GO TO SA 

59 UNITS = CUNITS 

NTI»1EP = T / 1AA0. ♦ IPNEXT 

IFIUNITS .FO. SFCS) MINER = T * 60. ♦ IPNEXT 
I FI UN I T S .EO. THIN) NTIMEP = T + IPMXT 
IFIUNITS .EO. HCUR1 NT IMEP = T / 6C * IPNEXT 
C 

NTIMFC= NTIMEP* NTIMEC - IPNEXT 
C 

C CCM/FRT FXC TO MATTS FCR GROOIN IF KiVENT ALREADY. 


1FUC0NVT 

.EO. 

0) 

G3 

TO 

650 

U«Z5 = C. 






IF t A I 315 I 

.LE. 

1.! 

I GO 

TO 

650 

D" f)05 JJ 

= 2, 

13 




JJ2 = J J 






IF(A(315| 

.LE. 

MATEXC 

I JJ 

,1)) GO TO R36 


8C5 CPMINUE 

UR 7 5 = 212. 

GO TO 650 

006 UR 2 5 = HATFXCIJJ2-1 ,21 ♦ ( ( ( A 1 3191-W AT EXC I J J2-1 , l > > / 

1 (HATEXCIJJ2, l) - WATEXCI JJ2-t,X) )', * IWATEXC(JJ2,2I - 

2 WATEXCI JJ2-X.2) )) 

C 

650 RETURN 

c 

C HERE IF DETFCTEO END OF RUN 18LANK TIME STEP CARO). 

66 UPZA s 3 
RFTU°N 
END 

SUBROUTINE FUNCTN(TH,PQL,TAB) 

0) ,/r NS ION TABI1A3 
N=1A 

CO 110 I=X,N,2 
IFUABUI-TH) 110,120,1X0 
110 CONTINUE 
GO TO 1A0 
120 PCL = T ABt I +1 } 

130 RETURN 
1 AO KN=N-2 





C'-i 150 I=l,NN, 2 

IPITARIU .LT. TH .AND. TABU+21 .GT . TH) GO 10 L60 
150 C J'lT I HUE 

WRITEI6.100) TH 

100 FOR M AT 1 5X t 8 ***** CURVE LIMITS EXCEEDED ***** ',G12.6//) 

IFITK ,LT. TAB! It > PCL«=TAB 1 2) 

IFITH ,GT. TABIN-ll) pcl=tabin) 

GO TO 130 

160 PQL=T AB t! +l)*(TABt 1*31— TAB (1+1 ))*lt TW-T AB(1))/{TAB(I+2I~TAB( DM 
GO TO 130 
END > 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE REMO <FUM , FUN2 , FUK3 , FUN4 I 
DIMENSION FUND l*»l .FUN2( 14) .FI1N3U4) , FIJN4 (14 I 

SUBROUTINE MEMO ( A«M ,4M3, ANU, ANY , ANZ , ARM , AUH , AUM , AUV, AVE , RFM ,BFN, 

* CN2.CN3.CN7.CV , GAS , DL A, DP A , DIU , DVS . F IS , HMD.HPL , 

* H*R .IISL ,'(S° . 12 . LVF PA ,PAM,Pa2,PC , PGL , PGS , PLA, 

* PPA.PP1.PP2, PRA.PR] ,PVS,QAO,OLN,OLO f OPO.QRF ,ORN, 

* OPO.OVO.RAM.RAR.KBF RPA , RPT , RP V , RSM, RSN , PVO.RVM , 

* RVS.U .VAE.VAS.VBO VI M, VI A , VIE , VP , VP A , VPC , VRA , 

* VPC .VPF.WE.VVR.VVS VV7.VV8.X , FUN1.FUN2.RUN3, 

* FUN4) 

REAL I2.LVM 

DIMENSION FUNIU4J , FUN2 114 1 , FUN 1 U4 ) . FUN4U4) 

REAL LVM,I,IFP,L'’U,KF,KE1,K0U,KIR,KI. ,KI , KCO, KEO, KN1 »KN3 
PEAL NAE.NFn.NID.NOO.n ,LPK,K tn.M0 2.IIOZ.KC7 .UPL.ttPR.I 2,I3.MMn 


COHMON/APRAY/T. I .VBO.VVS 

* VRE.PRA.'JRN 

* AIB.BPV.PPT 
COHHON/ARR AY/CN2 ,CN3 ,R V S 

* PAl.AUC.AUB 

* AU9,AUM,AU4 
COMMON /ARP AY /PC .Pf-O.VTC 

* 0PP.CHY.PPP 

* EPH.GP1.GP2 
CnHMON/ARR AY/R EK. NOD .NED 

* CNA.CCO.VID 

* OPZ.X .12 
COMMON/ ARRAY/ API, ARC, AP2 

* AN5.ANA.VB 

* PCD, VPC, RSN 
C0MM0N/ARRAY/AR2.PCC ,AP3 

* AHA.CNY.CNX 

* S° .VVk.RAK 
COMMON/ARRAY/LPK ,DPO ,HYL 

+ AHK.SRK.V9 

* HKR.PCV.PCZ 
common /arr a y/uhm . pod , 1 3 

* OSP.AftZ.AHY 

*• VPF.PPR.PPC 

COMMON /ARRAY/NOZ.KCZ ,VIZ 

* W9.02 A .01 

* PM0.P20 . PPO 
COMHON/ARR AY/Q 2 .03 ,PP5 

* PK3.FIS.STA 


VPA.VAS.VLA .VRA.VAG.PA.PAR.LVM, 

VPfc, PPA.PP1 .CPA.RPA.RV'*, VIE .PIA.OLN.PU . 
PGL, CPI', CS . WE. VVH ( PVS.PGV.R VG.OVO. AVE 
PGS.HTP.OAt , OR 0 , OLO ,DV S , CPA , CAS , PLA , 0° A, 

aun, aus.au;: ,A ua, oau.auj.au .auo.auh.vva, 
VIF.PPt ,PTT ,PTS ,PIF,C»I .PTC.CPP.PPC.PVG 
?t n.VTL.VTP.VPO.DPL.CPl.I.PC.nPI.t PI!, DIP, 
PGP.CPG.PCP.Gri.PGC.PGX.PGH.PG?. VGD.VG , 
GPO, AAR ,PK ,ftf N , APD.GL P . PFL . GFR , TRR , VU D 
NAE.VrC.CKj.KOD.RFl.KIO ,kIE,KCn,KE11 »CK[ . 
KE ,KI , VIC, II .VTY.Z .VIZ.VUZ.TVZ.PPZ, 
PP.l .VTS.VP ,PRP, ITP.GPR .KN3.KNI , AMR ,AMP 
AH3, AM5.AM ,CNE, AGK.ANP, AN1, ANC.AN2.AU3, 
HPl.HP ,Rl ,VIE,VI0,V|M,RC2 tPG 2, PKC.RCl, 
OVA.Uf N, 00. »,AOM,p 10,05 V, POT, PCD.POU, ARI 
ARM, CNI3.GF \ » AH7 , AHB, AH , AHC, AH( , AH2 , AH4, 
VVI ,VV2,VV5,VV6,W7,TVI)»VTH,HSR,HSL,N(G, 
CV .CNT.AUX.AUK.AUZ.Y .CFC.CPK.PCr ,cpr 
K in, AMT .AM .PUK.PtIN.AlK, A?K, A3K.CNB.CNZ, 
V2D.ZI ,Z2 ,23 ,Z4 ,25 , t-b ,27 ,ZR , HMK , 
ROO , CP2 , RBf , MC2 , POA , POY, ANU, PGR ,GF2 ,HtMO 
U , VP I , T 1 ,GF3 ,GF4 , AUP, AIJV, RV l , AUY .OUT , 
U5A , PP I , CP K , POS . PI F , PPG ,PPN, PRO , PF T , OF P , 
PMS, PMP.HR ,CPF,PfP,PAl,ilU, nOY.DPZ.GPZ 
W*P .HPL.ST, ,ALr,EXC»f2M,PA?,PP2,SVO,AUL, 
EXE , 4P F ,ORf ,PS ,JI ,BF u ,RA“,PVS,PVO,RMn.,onM, 
POO, PGE , API » A4K , PCM , CRM , P v 1 , PM 3. PM4 , CXI 
PK 1 ,29 ,Zli ,21 1,21.2, 213,21 4.2 15, Z16.PK2, 


PAR , GBL . ANT , ANZ , AIIX , ANV, ANm, ANfl ,AUO,AUR, 

* AUS.A378.HI ,AZ ,A3 ,01 MKIYU 9 I ,T( TLE 1 400 1 , DUMMY I40> 
COMHON/NUMFR(!/K,MG(20),NT I MFC , UN ITS. NZ.NT I ME P ,NN ,MAXNC, ATI ME 
C0MMGN/ST0RE/NGl,NG2 ,NG3 «NG4, NG5 ,NC< , NG7.NG8 ,NG9 ,OT,TLP,TNP,NO, 

* TM, TMM. NFIRST, 22(15), C.l YI9) ,0BYC9) , YMI N( 1 0) , YMAX ( 10 » , 

« N ,PT( IB > , BETA! 101 .NGRAI’HUO) .GRAPH (10 1 »HEARH91 ,NOEXP 



t 

JT 

U\ 

I . 



UTMAX 


C 

C CIMCKt ATfPY flYNAHICS fiLCCK 
C HfP'JQYNA'IICS 

c 

V9D = VO)-V'’ C-VVS-VAS-VLA-VPA-VRA 
VVS=VVS+CVS*I2*V3FI*.3986 
VPA=VPA+CPA*I2+vHn*. 155 
VA$=VAS+0AS*!2+VBD*.26l 
VLA=VLA+rLA*I2*VBO*.l28 
VP A=V" A(-DRA*r2*VBD*.C574 
VAF=V A S- .455 
PA=VA?/.0ih55 
■ IFfDA.LT.O.) PA-.0001 
P A u = 1 OG • /P A 
PA2=PA/AUP 

CALL FljNC TN(PA2,LVM,FUNl) 

VP P-V’ A- , 1 
PRA=V:>F/.005 

CAUL FUNCTN(PRA,QRN,FUN2) 

VPF=VPA-, 30625 

PPA=VPF./,C048 

PP1=.026*PPA 

!P|PP1.L7.0.) PP1»10.**1-12» 

PPA=PPlf*t-.S> 

PP2=°PA/ AUH 

[p| PP2.LE.0.) PP2-. 0001 
CALL FUNCTN (PP2 t RVM,FUN3> 

VL C=VLA- .4 
PLA=VLFF.01 

CALL FUNC1N (PL A.QLN. FUN4 1 

PPV = 1 ./IPLA+20. )/ * 0357 

RPT=RP V+R PA 

PGL='PPA-PLA 

CPO=PGL/PPT 

AUU=AT1 

IF (ANU.LT..8>ANU=.8 
VVF = VV«-VVP-(ANU-L.'»*4NY 
VV*=VVC-VV7 

IF ( VV6 .LT . .000 1 >VV8= .0001 

PVS=VV9/CV 

PR1 = PRA 

IF {PRA.IT.0.)PR1*=0. 

RVC,*2. 7RP/PVS 
cvo=(“vs-pan/pvG 

C';3=CN3M((PC-l7.)*CN7 + l7.)*CN2-CN3)v,l ‘ 
AVE-=(AU u .-l.}*AUY+l. 

RVS = AVF*U./CN3|*VIM+( (ANU-i.I*ANZ*l. I 
PG$=PA~PVS 

P SN=° AR A ARNPANUPAUMPPAM* VI M+R VS*1« 79 
RFN*PG5/F 5N 

ksm=anl*vifppam*aup#akp*hah 

PFNpPGS/PSN 

QAC=OFNFriFP+PPFKPA-PRA|*FIS 
OL'^LV'NAllLN^ AUH*H5L AHND+HPL 

Q?()*QBA* I t I •— QRF ) * AUH*RVH*HSR*HMD* HP H K}RF*QLO/QLN) 
QPn=OLCF l QPO-OLOJ /U 
QV3=tK0 + (CVG-0R0J/X 

dvs=oac-ovq 

opa«qro-qpo 


-9fi- 



* CAS =OLC-QAG 
DLA=OPC-CLC 
ORA=QVO-ORO 
RETURN 

EMI) 

SUBROUTINE AUTO IAU , AUB , AUC , AUH, AUJ ,/ UK, AUL »AUM , AUM, AUO, AUP ,AUQ , 

* AU D ,AUS,AUV,AUX,AUZ.IUA, AtJ6, AU8,AlB,OAUt EXC.EXE, 

* FX1.I2 ,PA t PAl,POO,fOT,P20,SrA,VVft,VV9,Y ,Z, 

* , Z8 ,Z12I 

REAL 12 

C 

C AUTCKC“tC CCNTRCL BLOCK 

c 

120 EXE = (B.-P20)*£XU(EXC-1.>*Z12 
POC=PCT 

If (PO<3.GT.B.)PnO = 8. 

IF IPCJ.LT.A.’)PCa=A. 

PAl=os*PCC/e.-EXE 

AUC=0. 

IFt PSt .LT.BO.).AUC=. 03*18 0. -PAD 
I F I ° \ 1 ,'t T .A0.)AUC = 1.2 
AU6-Q . 

IF(*A1.1T.170, |AUB=.0l*286*U70i-PAll 
IFtPAl.lT.AO. JAUB=1.83 
123 AlH=IAt)fl-l.l*AUX+l. 

12A AUN=0 

IFIP4t.tT.5C.)AUM = .2*(50.-PAl» 

IFIPAl.LT.20.)AUN=6.0 

AIJt = A t H-AU4 

AU? = A'JK<-( AU6-1 . I 

0AU=0ALMALCfAU6*AUN-DAUJ/Z/Y 

AUJ = MJJ + (nau-AUJI*I2*6./Z8 

IF(StlJ.LT.0.)AIJJ=O. 

IMAtJJ-l.J t26, 127,127 
126 AU=AUJ**AUZ 
GO TO 12 B 

12? AU=(*UJ-1.)*AUZ+1. 

123 IFISTA.CT ..0300UAU=STA 
AUC=AU-1. 

AUP=AU0«AU0*1. 

AIJH=AUC*AUV + 1 • 

AUR=AUC*AUS*1. 

VVR=VV9-AUL*AUP 

AUM=.15f.85*AUP 

RETURN 

END 

SUBROUTINE MORMON ( AM , AMC » AMP , AMR, AHT , AMI t ANM r CKE, P A, Z, FUN7, 

* AGK,ANC*ANP,ANR,ANT,ANV, ANW.ANl , CNA,CNE ,GFN, 

* I ,REK,A2 t A3) 

DIMENSION FIJN7UA) 

REAL I 
C 

c* * ** *o**#*4*» ****e***«********4*«**t*******4*"*************************»« 

c 

C ALDCSTERChF CCNTRCL BLOCK 
C 

c****^*******^**************************************************************^*** 

160 AMR-CKB/CNA/A3-A2 


I 

Ml 

I 



IF (m.LT.O. Mvci-o. 

CALI FUNCTN (PA.AMP, FLINT) 

AM 1= AM 1 + f ANM*AMP*AMR-4MU/Z 
AMC=A v CHAMl-APC)*(l.-EXP(-l/ AMT ) ) 

AM=20.Q39-I9.8»EXP{-.039l*AKC) 

C*******ft»»***#** 1.*4 *>!'*# 4 *:<t»**e ***4 >'*** 6 ** A********** ******** 

c 

C ANGIOT c NSIN CONTROL BLOCK 
C 

C ********■>**<■********* + ******************** ************************************* 

CNC=1 52.-C'Ja 
IFJCNF.LT. t. )CNE=l. 

AN»=I t 17.75-GFN*CNA)*AGK+1.)*REK 
ANVi=ANH* 1 IANR- t. ) * 10 .-ANW) AANVAJ 
1FIA‘M.LT..O)ANH=0. 

AN'P=ANR+ANW 

IFIANP.GT .100. JANPalOO. 

IFCANP.LT. .01) ANPs.Ol 

AN 1 = AN l ♦ ( ANP-AN1 )/Z 

ANC=ANC+( AM-ANC)*tl .-EXPI-I/ANT ) ) 

AN't=4.0-3.3»CXPC-.0967*ANC) 

IFIANP.LT..7)ANM=.7 

return 

END 

SUBROUTINE BLOOD <HKM,1-P ,HPK,I ,POT , P0Y,P01,P02,RC 1 .RC2.RCD ,RKC, 

* VB ,VIB,VIE.VIH,VP ,VRCI 
PEAL 1 

C 

C RED CELLS AND VISCOSITY BLOCK 

C BLCCD V I SCC S 1 TV 

170 VB=VP«VRC 

H1=100.*VRC/V8 
VIE*H W /(HMK-FW)/HKM 
VI e=V I F+l .5 
V I M= . 33 33*V 1 B 

C RED BLCCD CELLS 

PC2=PKC*VRC 

PI)?=PC1-PGT 

JF(P‘!?.LT..237 5)P02«.2375 

0C1=PCY*PC2 

PCD=»Cl-RC2 

VR C*VRC-* a CO*I 

RETURN 

END - 

SUBROUTINE PUSCLEC ALC , AMM, ACP, AUP, ARK , BFM, EXC, HM ,1 ,MMO,OMM.,OSAt 
A OVA,CVS,OZA,POO,PK1 , PK2, PK3 , PMO, PM1 , PH3 , PHA, PM5, 

* POE * PQH, PVOt P20i OO* 1 f RHO t VPF »Z5 Z6,RUZU 
RFAL I, PMC 

C 

C MUSCLF BLOOC FLOW CONTROL AND P02 BLOCK 
C 

180 GSA=AL0-VPF*.5 

C RUZ1- IS CXYGEN VOL. IN AORTIC BLOOOCCHBA) FROM GR0D1N. 

OVA « RUZ 1 * 1000. 



fifip non n o r> r> n o 


OVS^VS+t (BFH*CVfi-RHCI/HH/5./BFM-0VSJ;Z6 
PV0=5 1. 14*C!VS 

ru-is( PVO-PMOi «bms/(pmI**PK3-PM4 J 

W'=0CFMRPC-MMn)*<l.-EXP(-J/Z5)> 

PMj=PK2/(PKl-Qr]M) 

PMUPHQ 

IF(P M 1»LT «PN3 ) PNl =PM3 
P20=pmj 

IF(P2C.GT.8.)P2C=B. 

AOI*=(AUP-I.»#02A4l. 

KMO = Af:H*nH>1,'»EXC*n.-(8.0001-P20>P»3./3U.) 
poc=pvc-40. 

PkIE = PC'M*PQ04l, 

IF (PQE.LT ..005 IPO E=. 00 5 

A«m=AKM4 (PGE-AMM>*(1.-EXP(-I/A4K)1 

RETURN 

END 

SUBROUTINE AUTORG(AQP,ARK,AR1,AR2,AR3,AIK,A2K,A3K,BFN,DOB,HM , I, 

* M02,OSV,OVA,02M,PCA,POB,P(X,POD,POK,PON,PCR,POT, 

* PQV,P0Z,P10.0n2,PD0, Z ,Z4 ,Z7) 

REAL I,M02 

NCN-RIISCLF CXYGEN 0EL1VFRY BLOCK 

ANO NON-MUSCLE LOCAL BLUOO FLOW CONTROL BLOCK 


AUTnp.enuLATi cniRapio 


C$V=OSVH (PFN*GVA-0GBI/HM/5./BFN-0S\H/Z7 

POV=0SV*57. 1 A 

PDO=POT**3. 

IFCRQC.LT .50.)R0O=5O. 

01B-(POV-P(iTI«28S6,5/ROO 

MCi2 = Af>N*02M<-(l.-(8.000I-PtOI + *3./S12.J 

o<i2=or?+ (nrp-MC2i*u.-EXP(-i/ZAl} 

PUT-OQ2*. 00333 

PiO=PCT 

IF(D0T.GT.fl.)PI0=8. 

POC*PCV-P n P 

POa=PCBMPCK*PCD+l.-POBI/Z 

IF(POB.LT..2)P0R=.2 

ARl = ii>lt(P0D-6Rt )*( 1.-EXPI-I/A1K) I 

ARH = A* > l*A l5 ?*A c 3 


AUTORtGULATTCN, IM FP"FD I ATE 

pnA=POA+(PON«-POO + l.-POA]/Z 

IF (PCA.LT..5IPCA=.5 

AP 2=A R 2+ l POA— AR2 J*ll.-EXP(— I/A2K1) 


AUTORFGUL AT I CK , LCNG-T ERM 


IFIPOSI194, 192,192 
192 P0C=PC2*PC0+1. 

GO TO 196 

194 P0C=PCZ*PCC*.334l. 

196 IF [PUC.L7..3)P0C».3 

AR 3=AR3* ( P0C-AR3 > + 1/ ASK 

RETURN 

END 


1 

NJ1 

VO 

1 



SUBROUTINE ADH {AH , AKC, AHK , AHM, AHY , AHZ, AH7, AH8 , AUP, CNA.CNB, CNR» 

* CNZ,t ,PRA,Z) 

■ PEAL I 
C 

C ANTIOIURETIC HORMONE 

c 

CNB=CNA-CNR 
AUZ= . 2*PR A 

AHY=AHY+( At-Z-AHY)*.0007*I 
AHB = AUP-1.. 

IF(AHP.LT.0.)AH8=0. 

1F(CN3.LT.0.)CNB=0. 

AH=AK+{ CNZ*CN8+AHB-AHZ+AHY”AH)/ Z 
IF(AH.LT.O.)AH=0. 

AHC=At-CM •3333*AH-AHC) * ( l.-EXP (-I/AHK ( ) 

AHM = 6.4( l.-EXP{-0. 1808*AHC)) 

IF{AHP.LT..3>AHK=,3 

RETURN 

END 

SUBROUTINE HISC1 { AHM, AUA , AU , I ,SR v !>RK , STH, T VO, TVZ, VEC . V 1C , VTW, 

+ VVE f VV 6 1 VV,7 ,Z tV9' 

REAL 1 
C 

C*44* ****** 4 ************** 44»****4***4*44*:<'>i4*444*44****4**«*44*44*******#4*4*** 
C 

C VASCULAR STRESS RELAXATION BLOCK 
C 

£*4*444*44 4 444 4 4* *444 *444* *44* 44**4* 44 4*4*<.'.»444***44****44 4 4*4**44*4 44 4***a*4*ff 

VV6=VV£>AtS®*(VVF-V9 1-VV7-VV6J/Z 
VV7=VV7<VV6*( l.-EXPC-l/SRK I ) 

£444*444**44** 4 4 4 444 4 44444 44 4444444*4 4444 4* •<444*44*444444*44 44*4*4*44 444*4484*4* 

c 

C THIRST AND CR INKING BLOCK 

C» 444**4 44 4 4 44 444444**4*9 **»***a**4**44***»- t********* **44*44*4*****»*4 ********** 

c 

TVZ = { .9 l * AH W - .009 ) *S T H 
TVO=Tvn*(TVZ-TVO)/Z 
IF|TVO.LT.O.)TVD*C. 
vn.=vtc+VEC 

C 4*4* ****** 4 4*4 44**444**4* *4**44444 44*44*4 4 (4 4444*4444*44 44**4*44*444*4 444**4*44 

c 

C AUTCNCWJC CCNTRCl PLCCK 

C 40APTA 1 ION OF BAR0RECEP70RS 

C 

c* *44*4 * **4 44***444***444*4* *44*4444*494444'I4444*444444***44444**444**444*4*444* 

AUA=AU4*AUB*T 

RETURN 

EN9 


C 

c 

c 

c- 

c 

c- 


* 


SUBROUTINE HEART IAUR .CHM.HMD.HR .1 .PA- *PMS rPMP.PMSaPETaPRArO'Afl* 
OLE .R7P»SVO.VAE»V6.L .VPE ,VRE rWEi-KIX 


REAL i 


HEART HYPERTROPHY OR EFTERtORAT IGN- BLOCK 


HEART YSCI'CUS CYCLE 



OH**! PCT-6. )*.0025 
HMR = H u, )+OHMAI 
IF H'G.GT.U}H , 'n=l. 

C MEAN CIRCULATORY PRESSURES 

PMC=( V»E+VVE+VAE«-VPE*VLE)/.ll 
PMS=IVAE + VVE+VRE)/.0*n75 
P«P=< VPEtVLE)/. 0162b 

C** = *»*» *■»»#= = = = # = *#=»«»=**** = »**»**« «*»»■»*« »*»*=*»** o =v = »****r*=« COUO «**«=**« O* 

C < 

C HEART RATE AND STROKE VOLUME BLOCK AND TOTAL PERIPHERAL RESISTANCE 
C 

C**»A ******* ******$***»»#*$♦*'>«*.♦*'<':*'!'**♦*<'#*« ♦♦**++****♦**** #*****#****#«. #**(■#** * 
HR = t 32.+H1 «AUR+PRA*2. I'M t HMO-1 . I*.5+ 1 . > 

PTP = ( P A— PR A ) /Q AQ 

svc=olg/hr 

RETURN 

ENO 

SUBROUTINE CAPMBDC BFNrCFC.CPI ,CPP,DFP, l ,IFP»PC ,PCD t P IF,PLO ,PPC * 

* PRP« PTCtPIS t PTT » PVGi PVStRVStTVDf VG ,VIO,VIF,VP* 

* VPO,VTC,VTO,VTL t VTS. VU0,2 ,Zl,FUN6> 

DIMENSION FUN6I1A) 

REAL I , IFP 
C 

C CADHLA»Y MEMBRANE DYNAMICS BLOCK 
C 

130 PTTMVTS/12.)*»Z. 

VIF=VT S-VG 

CALL FUHCTN IVIF f PTS,FUN6l 

PIF=PTT-PTS 

CP I = t FR/V IF 

PTC=.25*CPI 

CPP=PRP/VP 

PPC= .AACPR 

PVGsOVS*! .7S#BFN 

PC=PVG*PVS 

PCD=PC+PTC-PPC-PIF 

VTC=VTCMCFC*PCD-VTCJ/Z 

PL 0=7 .tHPIF— PTT 

VTL=VTL+( .OOA^PLO— VTL 7/Z 

IMVU.LT.O.) VTL=0. 

VTD=VTC-VTL”VID 

VTS=VTS*VTD*I 

VPD=VPD-f i TVO- VTC+VTL— VUU~DFP—VPD)/Z1 

RETURN 

END 

SUBROUTINE PULMONI CPF*CPP>CPN< DF'P*l <'PCP» l PFH^PtA>PtFYPOSYPPA»PPCY 

* PPOyPPH'P'f'Nv-P'PCv’P’PP* VP »VP'0 >VP'Fv4 >234 
RPAL 1 

c 

C PULMONARY DYNAMICS AND •F‘LMjt‘05 BL'PCK 

c 

VP=VP*rtVP0*T1/Z3 

■c 

200 PCP=» A5*PPA+«55+PLA 
PP1=2.-„<IS0/VPF 
'CPN=PPR / /VPF 



PPS=CPA*.4 

PLF=(PPI+I1.I*.0003 

ppri=PLF*CPN 

P°N= ( CPP-CPN }* «OOQ225 

PPf) = PP0* < PPN-PPn-PPO ) U 

IF(PPRtPPl)*I— .C2 5»LT.0«JPPO=( .025— P PR J /I 

PF I = IPCP-PPIfPOS-PPC)*CPF 

DFP=QFP»( PFI-PLF-DFP) /Z 

IF (VPF»CFP*I-.001.LT.O.IDFP=( .OOl-VPFJ /| 

VPF=VPF+t)FP*| 

PPP.=PPR*PP0*I 
Pr ri/' 1 

tuu 

SUBROUTINE MISC2 CHPL,HPR,H$L ,HSR, J,Pt .PPA ,POT, STH.ZlO* Zl 1 .21 3) 

REAL I 
C 

£*********************»***<****************<. < 4*9 ************ 4 ************* ****** 

C 

C PEART HYPERTROPHY cr ceterioration block 

c 

C**»* ***«**»»* ******* ************** ********* >*********************************** 
HPL = HPL*t ((RA/IOO./HSL )**Z 1 3) -HPL}*I / 3 T600. 

HPP=HPR*( ( tPPA/15./HSR)**ZL3l-HPR)*I/ST60C. 

C ♦ *** *** *** * ** * ******************** ********* ******************* A*** ************* 

c 

C TISSUE EFFECT ON THIRST AND SALT INTAKE 
C 

£************** ********** ******************* >*********************************** 
STH=( ZIO-PCT) *211 
IFCSTH.LT.l.)STH=l. 

IF(STH.GT.8.1STH=B. 

RETURN 

ENO 

SUBROUTINE PROTFN<CHY,CPG,CPI,CPK,CPP.CPR,CPl,OLP»OLi,CPC.DPI,DPLt 

* UPC* CPYtGPO »GPR* I ,IFP,LPK,PC ,PC6,PGX,PRP,VG . 

* VTl.Z ,PPD) 

REAL I > I FP »LPK 

C 

C TISSUF FLU IPS iPRESSUR FS AND GEL BLCCK 
C 

C PLASMA A NO TISSUE FLUIO PROTEIN 

135 OPL*OPLitVTL*CPI-DPL 1/2 
IF CPC.LT.O.IPC=0. 

COC^ilRC+lCPKR-ICPP-CPI )*PC**PCE-OPCJ/Z 

OP l =QPC-G 3 L 

DLZ»LPS* (CPR-CPP1 

IF tCPP.GT .CPK )f)L2*A.«0LZ 

HLP=OLR*t C1.2-0LP1/Z 

pap=PRPf ICLP-OPO+OPL-CPC-PPOI*! 

C . GEL PRCTEIN GYNAVjcs 

141 PGX*=CHY**2*.01332*CPG+CPG 

GPD=GPO«- ( . 0005* CCPI-PGXI *VG-GPD t/Z 

GPR=GPK*GPD*I 

IFP«1FP+(CPI-GPD)*I 



non 


nf turn 

END 

SUBROUTINE KIDNEY! AAR, AHM, AM , APO, ARr , AIM,CN£ ,CNX,CNY , GBL , CFN, GFR , 

* GF2,GF3,GF4,GLP,I ,NAE , NE0',NI 0,N0D, NOZ.PA ,PAR, 

* PFL.PPC,RBF,REK,RFS,RR ,STH t TRR,VIM,VUD,ZI 
REAL I ,NAF , NED ,NI0 ,NG0 ,NOZ 

c 

C KIONFY DYNAMICS AND EXCRETION SLOCK 
C 

162 GF3=( (GFN/. 125-1. I*GF4)*1. 

IF IGF3.GT.IS.)GF3=15. 

IFIGF3.LT. • 41 GF3=«6 

AAR = 31 .67*VIH*( AIJM*ARF+1.-ARFI*GF3 

PR=AAR + 51 .66*VIH 

PAFbPA-GGL 

PFN=PAR/PR 
R RF = R EK*P FN 
150 APD'AAR’tP FN 
GLF'PAR-APD 
PFL=GLP-PPC-18. 

GFl=GFN 

GF N=G FM t PFL«.OQTPI-GFN) *GF2/Z 
IF (AaS(GFN~GFn.GT..002»GO TO 152 
GFRsGPN*KEK 

TPR=. B-»GFP + .02 5*RFK-.001»RFK/AM/AHH 
VUC=VUD+ ( GFft-TRR-VUO >/Z 
IFIYUC.LT ..0002 >VUO=.0OO2 

C K1DNFY SALT OUTPUT ANC SALT INTAKE 
C I SEE ALSC ELECTROLYTES ANO CELL WATER BLOCK ) 

noz=iooo.*vuo/ak/ccne/cnx+cny> 

M0D=NI1DyI NOZ-NOO) /Z 

NEC=NIC*STH-N0D 

NAF»NAE«NED*I 

RETURN 

END 

SUBROUTINE IONS I AM , CCC.CKE ,CK I ,CH, I ,KCO,KE ,KE0,K1 ,KID,K1E, 

* KlR,KOO,NAE,REK,Vt:,VIC,VID,VP ,VPF,VTS,Z> 

PEAL I ,KCD,KE.KEO,KI ,KID,KIE*KIR,KOC,NAE 

ELFCTROLYTFS ANC CELL HATER BLOCK 

160 VFC=VTS«VP+VPF 
CKFsKE/VEC 

KCD=l . 00042*CKE*.00014*AH*CKt»*REK 
KIR = 2R50. *UO.*CKE 
K I F-K IP-K I 

KCD=KCO+ I K I F». 013-KCC) /Z 
K !=K t+KCn*I 

kec=k.ic-kco-kqd 

KE*KE tKEO *1 
CKI=K I/VIC 
CNA»N4E/VEC 
CCO=C K I-CNA 

VIC=VIO-*( .01*CCD-VI01/Z 

VIC=VIC*VIO*I 

RETURN 


I 

CT\ 

(JO 

1 
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• END 


SWO®TJTINE GELFLD(CHY,CPG,CPI,GPR,HYl . 1 FP , PGC , PGH , PGP, PGR, PGX, PIf , 
* PRM, PTC,PTS.PTT,VC « VGO, VI F , VRS, VTSt VZD. FUN6) 

OIHFNSION FUN6I14) 

REAL I FP 

GEL FLUID DYNAMICS 
140 CHY^HYL/VG 

PRV=-5,<j*CHY+24.Z 

PGR -• 4^00 Y • 

CPG'GPO/VG 
PGP=.2B*PGX 
PGC=PGPFP5P 
V I F = T S-VG 

CALL PUKCTN ( V I F tPTS «FUN6I 

P|F=PTT-PTS 

CP[=IFP/VIF 

PTC = .2S*CP! 

PGH=PIF+PTS+PRM 
VGO=VZO*IPIF f PGC— PTC-PGH I 
VG=VG«VGD 
IFIVG.LT. 0.)VG=0. 

|F{.012.LT.A8S(VG0I ) GC TO 140 

RETURN 

END 

SUBROUTINE GPOOINIURZI,URZ2,URZ3,URZ4,URZ5,RUZ1I 


FOLLOWING FROM GUYTON TO GR001N 


GUYTCN 

GRCOIN 



OLC 

CUO) 

(URZ1 ) 

CARDIAC OUTPUT. 

(RMO+DUB ) /10Q0. 

RWT l 2 )+CI 26) 

(URZ2) 

TOTAL METABOLIC RATE OF BODY 

(HM*OSAI/200. 

C 1 17) 

(URZ3) 

BLOOD OXYGEN CAPACITY. 

URZ4 

URZ4 

CURZ4I 

FLG 1=1ST.TI*E GROOIN CALLED 




2=>GUYT0N JUST OUTPUT 
3 s GUY ION JUST OUTPUT 
ANO DETECTED END RUN. 
0=N0T ABOVE. 

URZ5 

URZ5 

IURZ5 ) 

WORK LEVELIWATTS ) . 


FOLLOWING FRCM GRCOIN TC GUYTON 

GUYTON GROOIN 

OVA/ 1000 . CHBA (RUZ1) OXYGEN VOL. IN AORTIC BLOOD. 

Cl 40 ) 

ALVEOLAR VCL GAS FUNCTIONS 

1 FAICC2I 

2 F A I 07) 

3 FAIN2) 

GAS CONCENTRATIONS IN BRAIN. 

4 CBICC2) 

5 01(02) 

6 CBIN2) 

GAS CONCENTRATIONS IN TISSUE. 

7 CTIC02) 

B CTI02) 


~*i9- 



c 9 ctinpi 

C CARDIAC OUTPUT. 

C 10 0 

c cerebral aiono flow. 
c u aa 

C GAS TENSION IN CSF. 

C 12 PCSFICC2) 

C 13 PCSFIC2) 

C 14 PCSFIN21 

^ . 

C LFNGTH OF SIMULATION RUN. 

C 15 TMAX 

C WEIGHTING OF H+CONC IN CSF VERSUS VENOUS BLOOD OF BRAIN. 

C 1ft CENTRAL SENSITIVITY PARTITION 

C BL COO OXYGEN CAPACITY 
C 17 (H0) 

C time CONSTANTS IN CARDIAC OUTPUT AND CEREBRAL BLOOD FLOW RESPONSES. 
C IB Rl 

C 19 «2 

C 

C CONTROLLER ECUATIGN SENSITIVITY WEIGHTINGS. 

C 20 CENTRAL SENSITIVITY COEFFICIENT 

C 21 CAROTID BOGY SENSITIVITY COEFFICIENT 

C 

C VOLUMES CF LUNG » BRA I N , AND TISSUE 
C 22 KL 

C 23 *9 

C 24 KT 

C 

C BRAIN METABOLIC RATE OF C02 PRODUCTION. 

C 25 HRBIC02) 

C BRAIN METABOLIC RATE OF 02 CONSUMPTION • 

C 26 RRBIC2) 

C GAS DIFFUSION COEFF.FOR BLOOD-BRAIN -BARF IER. 

C 27 0CG2 

C 29 002 

C 29 CN2 

C 

C BAROMETRIC PRFSSURE. 

C 30 1! 

C VOL. FRACTION CF INSPIRE!) GAS. 

C 31 FIIC021 

C 32 F 1 102 1 

C 33 FIIN2) 

C 

C VOL .OF CSF. 

C 34 KCSF 

C INITIAL TIHF 

C 35 T 

C COMPUTER TIMC STEP. 

C 3ft H 

C CONTROLLER FDUATIOK CONST ANTIHAINTAINS NESTING PAIC02I APPROX. 40). 

C 37 VI I N ) 

C VALUE FOR RESTING ALVEOLAR VENTILATION. 

C 3B VI (55) 

C OUTPUT PRINT INCREMENTS (ALSO PRINTS AT .5MIN. INCRIMENTS) . 

C 39 PRINT-ALL TIME 

C 

C SV t 18.50) 

C ARTERIAL GAS CCNCENTRAT IONS AT LUNG EXIT. 


-S9- 
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1 CA(C02> 

2 CAl 02 > 

3 CA(N2) 

VENOUS GAS CCNCFNTBAT IflNS AT BRAIN EXIT. 

A CVR(C02) 

5 CVEIP2) 

6 CV6IN2) 

VENOUS GAS CONCENTRATIONS AT TISSUE EXIT. 

7 CVTIC02) 

B CVTI02I 
9 CVTIN2) 

CARDIAC OUTPUT. 

10 C 

CEREBRAL BLCCD FLCH. 

11 OB 

TISSUE BLCCO FLOW. 

12 CT 

ARTERIAL hi CONCENTRATION. 

13 CA(H+) 

ARTERIAL 02 TENSION. 

I A PA ( 02 J 

15 — 

TOTAL GAS CONCENTRATIONS AT BRAIN EXIT. 

lb CVBIC02) 1 CVBI C2 I + CVBIN2) 

TOTAL GAS CONCENTRATIONS AT TISSUE EXIT. 

17 CVTIC02J ♦ CVTIG2) + CVTIN2J 
TI^F. 

18 T 

' VTRANUB) 

ARTERIAL GAS CONCENTRATIONS AT BRAIN ENTRANCE. 

1 CAu [ G02 1 a CA(C02H T — TAB) 

2 CABt C2 ) * CA102IIT - TAB) 

3 CAR ( N2 ) = CA I N2) I T - TAB) 

VFNCUS RRAIN GAS CONCENTRATION AT LUNG ENTRANCE. 

A CVR (CO?) { T— TVB) 

5 CVO (C2 ) ( T - TVB ) 

6 CVBIN?) (T - TVB) 

VENOUS TISSUE GAS CONCENTRATION AT LUNG INTRANCE. 

7 CVT ( C02) ( T~ T VT ) 

B CVT (02) ( T - TVT) 

9 CVT(N2)(T - TVT) 

ARTFPIAL GAS CONCENTRATIONS AT TISSUE EN RANCE. 

10 CATIC02I = CA(CC2)(T - TAT) 

11 CAT 1 02 1 = CA 1 02 ) ( T - TAT) 

12 CATCN2I = CA(N2)(T - TAT) 

ARTERIAL H+ C0NC6NTRATICN AT CAROTID BOU ES'SITE. 

13 CAO(Hl) = C A ( HI ) ( T - T40) 

ARTERIAL 02 TENSION AT CAROTID BODIES'S)' E. 

19 RAO (02 ) - PA(021(T - TAO) 

ARTERIAL H+ CONCENTRATION AT BRAIN ENTRANCE. 

15 CAft(Hi) * CAtH+MT - TAB) 

TOTAL GAS CONCENTRATION FROM BRAIN AT LING ENTRANCE. 


t 

On 

on 

i 
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16 < CVB tCt 2) * CVM02) ♦ CVEINZ) ) U-'TVfl ) 

TOT 41 A 5 CONCEMTRAT ION FROM TISSUE AT UING ENTRANCE. 

17 (CVT(Cr2» ♦ CVTI02) + CVT(N2))U - TVTJ 

CI15 J 

FOR Cl 1 5 > The SYMBOLS 0=8ARCMETRIC PRESl.URE, 47=WATER vapor press., 
K=CONVERSIQN FACTOR FOR ATM TO MMHG * A=5<1LUBI LI TY CCEFF.OF GASES, 
H=COMPUTFR T|«E STEP, HB=BLCOO OXYGEN CAPACITY 

1 B - 47 

2 K AC02 

3 K A02* 

A K 4N2 

5 K AN 2 IB - AT) 

6 K A 02 IB - 47) 

7 K AN Z <B - 47) 

8 0.16 t 2.3IH6I 

9 863/ ( P - 4 7) 

10 0.62 

11 K ACSFIC02) 

12 K ACSFI02) 

13 K ACSr|N2l 

14 2*M 

15 1. 99*H 
F I 20) 

compartmental gas tensions and concentra'ions. 

1 PAI02) 

2 K AC 32 PAIC02) 

3 P3IC2) 

4 K AC02 PBIC02) 

5 PTI02) 

6 K ACC2 PTIC02) 

7 PAIC02) 

8 PA 1 02) 

9 C A { 0 2 ) 

10 CA tti2 ) 

11 CAIC02) + CA (02 I ♦ CA(NZ) 

12 CVBI021 

13 CVTIC?) 

PRODUCT OF O | FFUSI CN COEFFS.ANO GAS DIFFERENTIALS ACROSS BLOOD-GRAIN 
BARRIER. 

14 QCC2 (P8ICC2) - PCSF( C02 ) 1 

15 D02 |PP<02) - PC5F10Z) ) 

16 DA2 (PDIN2) - PCSF( N2 ) ) 

17 PE (02 ) 

15 Pfl(N?) 

INTCCER UR24 

DIMENSION C (40 1 , XN(40,2>, SVU8.50), VTRANU8), RKI14.4), 

1 SC( 14,5) , CC(I4), A[6 ) , 0(15), F(20), VOL(IO), RMT (2) , 

2 BC14), 0 F ( 6 ) , TAU ( 5 ) , CCl3i, CHB(3), CH(4), CPMI3), 

3 D0(41 

DIMENSION XNBI4.2), 0JC4), C0J(6), I0JI2I 

COMMON/2/ C, XN, SV, VTRAN, RK, SC, OS, A, D, F, VOL, RMT, BCi OP, 

1 TAU, CC, CUB , CM, CPU, DQ, VE, VI, CPB, CPT, C ADK, X, DT, 

2 IRK, LOC, ITERX, INDEX, 1, J, 4, N 
CUMMON/P/ XDS,XMH,CXT, WORK, DUMl, 0UM2, 3UM3 ,H0RK2 ,RMTB, RMTB2 ,T IMEOF 

1 ,RHLIN 
C 

c 


» 

o\ 


1 



C SEF IF THIS 1ST TIME GRODIN CALLED. 

IMUR2A • NE« I I GC TO 370 
C 

c 

C HERE IF THIS 1ST T IMF GRODIN CALLED. 

C 

C DATA FnR INITIAL CONDITIONS 

WRITE (6,5) 

5 FORMAT (/» GRCDINS: RFSP I RATOR Y COMRCL MODEL 1 ) 

300 C n hT INUE 

WPITE(6,90) 

SO FORMAT ( 1 0 RESPIRATORY CHEMOSTAT — If PUT DATA'/) 

C DATA FOP INITIAL CONDITIONS 
DH 10 I = 1 , AO 

RFAD(5,19C,FND=30I) C t I ) , I XN t I , JJ ,J=l,2l 
10 CONTINUE 

c establish, computer step independent of input data. 

C FOLLOWING TIME STEP TO MAKE GRODIN COMPATABLE TO GUYTON. 

CI3&) = .003 

ISO FORMAT (.SX,F15.0,5X,2AA) 

CO 20 I = I, A 
IPAD a I ♦ AO 

READ (5,190) BC ( 1 1 , (XNftU.J), J « 1,2) 

20 CONTINUE 

DO 30 I = 1,2 

PFAD (5,190) RMTl I) , ( XNB( I , J) , J 3 1,2) 

IP AO = I t H 
30 CONTI NUF 1 
DU AO 1 = 1,2 

RCAIJ 15,1901 DJ ( I ) , (XNB(IiJ), J - 1,2) 

IPAO = I ♦ A6 
AO CONTINUE 
C 

c 

CHOI * UR21 
C ( 1 7) 3 UP23 
RMTl 2) « UP22 - Ct 26) 

C OUTPUT INPUT ♦ ABOVF 3 VARIABLES FROM GLYTON. 

J 3 1 

00 75 I - 1,8 
JJ 3 J+A 

HR I TE I 6 , 92) J, (CIII).I1«J,JJI 
92 FORMAT,! • ' , 12 , 2X, 5t F9 . A ) ) 

J « J ♦ 5 
75 CONTINUE 
WRITF(6,92) 

J » A5 
H'f I TF ( 6,92) 

C 

C FKC02) 

DUN lsC( 31 ) 

C FKC2) 

01)92*01 32) 

C FI (7)2 ) 

DUM3=C(33 > 

W'TRKxO, 

H0PK2=0. 

C METABOLIC RATE 
RM T8»RMT (2) 


J, (Bcm,!-X,A) 

J, RMTl 1 1 , RMT(2), OJ(l), DJ (2 ) 


OF C2 CONSUMPTION IN TISSLE. 



no no n non 


RHTU2ssRMT(2> 

C 

T [‘•■POP -C • 

X»S=0. 

X*H = tO. * CO6I/.0C3 
H'tMsrO 

201 CONTINUE 

X>3S=X0S<-X W H 

I F ( HMH .£0 . t )XOS=XOS+C( 36 ) 

NMM = 1 

C t 35) = C« 

cuoi>o. 

INITIAL GUESSES FOR ITERATIVE LOOPS 
ARTERIAL CONCENTRATION OF C02. 

CC(1) = 0,6 

BRAIN CONCENTRATION OF CC2. 

CC 12 ) = CIA) 

BRAIN C 32 TENSION. 

CPS => 50.0 
TISSUE C02 TENSION. 

CP T = SO.C 

IFIXDS.GT -XMH) GCT0202 

SFTS VARIOUS CONSTANTS AND AGGREGATES OF CONSTANTS- 
TMAX. 

CHS) *= CU5) ♦ .0001 
C PRINT ALL TINE. 

CIST) = C(39) + .0C01 
DO 200 I - 27,29 

C FACTOR OF l-E-7 MULTIPLYING DIFFUSION COt FFIC IENTS. 

cm = cm * i.E-T 

200 CONTINUE 
202 CONTINUE 
I?K * 1 
Ms* 1 A 
N s* 5 

lujtn = c 

C SCIU-JILI IT COEFFICIENTS. 

C A(I)= (ALPMAJC02, A(2> = (ALPHAIQ2, A(3)» (ALPHAJN2, 

C A( A ) = ULPH4)C02, A (5 I *= t ALPHAI02, At6l * IALPHAIN2 

■ AIII = 0.51 
AI2) = 0.02A 
A ( 3) » 0.013 
A I A ) = C. 5 1 
4(51 = 0.92A 
At 6) = 0.C13 

c AT M/m -I HG CONVERSION FACTOR. 

SK = 0.00132 

C CARBONIC ACIO nlSSCClATIDN CONSTANT. 

CACK = 795.0 

C VOL 1 1 1— VOL 1101= VCLUHES USED IN CALCULATION OF VARIABLE TIME DELAYS. 
VOL ( 1 ) = 0.015 
VCLC2J = I.C62 
VJLI3) = 0.1 93 
VOL I A I * 0.06 
VOL I 5 ) * 0.109 
VOL ( 6 1 = 2.9A 
VOL I 7 1 *t 0.735 
VOL 18) = 1.062 
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VOUJ) 3 0.003 
VULU0I = 1.062 

MET ABOL 1C RATE OF C02 IN BRAIN ♦ TISSUE. I / SAME FOR 02 
QR<6) = (C(2S) ♦ HHT(l)J/(C(26) ♦ RHI12IJ 
a-47 

0 1 1 1 = C ( 30 I— 47 • 

210 I = 2,4 

PRODUCTS CF COJVERSIflN FACTORS AND SOLUBILITY COEFFICIENTS, 
nil) = SK*A(t-U 
DIH9I a SH*A I | +2 ) 

0(1*3) = D(l)*DU) 

210 CONTINUE 

FACTOR USED IN ESTABLISHING CAICH2I 
0(8) = 0.16 + 2.3KU7J 

0(9) = 063. 0/0(1) 

FACTOR USE0 IN ESTABLISHING CBIC02). 

0(10) = 0.62 

KAMPULAT ICN CF COMPUTER TIME STEP. 

0(14) = C ( 3 6) *2. 0 

0(15) = 0(141 - .01*0(36) 

CALL RC3 
CALL RC4 

call KC5 (CPR, F(4>, C(4), BCI21) 

CALL RC21 ( CH B ( 2 ) , F(3), F(4), C(4), (H(2), CPH(2)) 

CALL KC 19 (CPB, CHB( 2) , CC(2), BC(l), Ft4>) 

CALL RC5 (CPT, F (6 ) i C17I, 6C<3)) 

CAIL RC21 ( CHB ( 3) , FIS), F{6>, C<7), (.H(3), CPH(3)) 

CALL RC19 (CPT, CHB(3), CC(3), BC(1), F(6)> 

CALL MC20 
CALL RC7 
CALL RCB 
CALL RC9 
CALL RC10 
CALL PCI l 
FORCE PRINTOUT 
URZ4 = 2. 

CALL RC12(URZ4,UR25) 

GO TO 60 
50 CALL PCI 5 
CALL RC 1 6 


60 RUZ I = CH 0 1 1 ) 
RETURN 


NORMAL ENTRY FROH GUYTON. 

370 ctio) 3 urn 

C( 17) = URZ3 

RMT ( 2 ) = URZ2 - CI26) 

CALL KCl 3 

CALL ftCl2(l/RZ4,URZ5l 


SEE. IF. GUYTON SAYS END OF RUN. 
IF1URZ4 .64. 3 ) GO TO 6G 


I 

■'7 

O 


I 



IF 1C I 35) »G8 .XMH) GC TO 201 
C 

73 CALL PC14 

UU = AMPO(C1351. 0(1411 

IF CUU .LT. .0001 .OR. UU .GT. 0(15)1 ' G0T050 

GOTO 60 
C 

C HERE WHEN GUYTON SAID END OF RUN. 

80 IF (Cl 37 ) .GT. 1.0E-5) GO TO 250 
220 CTFRM = 0.0 

IF (VTRAN ( 14 ) - 104.0 230, 240, ,40 

230 CTER’1 = (23.6E-9) A( ( 104.0 - VTR AN ( 14) 1 **4. 9) 

240 CC 371 = C(20)*1C( 16) *VTRAN(15) * (1.0 - C ( 16) ) *CH(4I ) 
l ♦ C(21)*VTRAN( 13) ♦ CTERM - VI 

I = 3 7 

WRITE (6, 192)1, C( I), (XNU.J), J = 1,2] 

250 DO 260 I = 1,14 

W<ITE(6, 192)1, C(I), (XNtl.J), J * 1,2} 

260 CON T INUt 
301 CONTI NUF 
STOP 

192 FORMAT ( 13, 2XF15.5,5X2A8) 

194 FORMAT (lHl) 

FND 

SUBROUTINE RCB 

DIMENSION C (40 ) , XN(40,2), SV(13,50), VTR AN 1 10) , RMI4.4), 

1 SCI 14,5 ) , OC (14), A (6 ) , Dili), F(20), VIILUO), RMT12) , 

2 RCI 4) , OF { 6 ) , TAU( 5) , CCI3), CH0I3), CH14), CPH(3), 

3 0014) 

COMH0N/Z/ C, XN, SV, VTRAN, RK, SC, CC , A, D, F, VOL, RMT , 0C, OF, 

1 TAU, CC, CH6 , CH, CPN , 00, VE, VI, CPB, CPT, CAOK, X, OT, 

2 IPK, LOC, ITFRX, INDEX, I, 4, *, N 

C CALCULATES TRANSPORT TIMES 

C ECUATIUNS B.IO THRU S.14 . 

C6969 F0RMAT11H 7HSUB RCB) 

00 B70 I = 1,5 

CT = C( 35) - SVtla.l) 

NO = 1 

GO TO (810, 812,814,816,810), I 


310 

NC 

= 11 




NR 

= 1C 




r,n 

TO 820 



012 

NC 

= IC 




NO 

= 11 




GO 

TO 820 



814 

NC 

= 10 




MR 

= 12 




GO 

TO 820 



816 

NC 

= 12 




NR 

= 10 




CA 

= OF(l) 




GO 

TO 822 



820 

QA 

* C1NC1 



822 

00 

860 J * 1,2 




GO 

TO (834,824), 

J 


624 

NC 

- N8 




NO 

« K + 1 




IF 

IK) 

826, 826, 

832 

826 

IF 

INC - 12) 

830, 828, 

830 


1 

^) 

I 



OUJU 


62* 0 4 = SVINC.U - ($V(NC,l) - OF 1 1 ) )*CT ' (C( 35} - SVU8»1I) 

00 TO 834 

830 OA = SV(NC,L> - ($V(NC,I> - C(NC) )*DT '(C(35) - SVU8, 1)) 

GO TO 83 A 

832 OA = SVINC.ND) - (SVtNC.K) - SV(NC,NU) )4QT/D(14> 

834 IJ * 2*1 ♦ J - 2 

Ad = VOltlJ) * (C(36)/. 0078125} 

AA s 0T e ( OA ♦ SV(NC.NO) >/2.0 
DO 833 K = NO, 49 
IF (AA - AB,1 ' 836, 836, 840 
036 AA = AA «• C(36)*(SV(NC,K) * SV(NC,K+1)J 
838 COST INUE 
K *> 49 

WRITE <6,e90I I 
840 DA = AA- AB 
K = K -1 

IF (K) 842, 842, 846 

842 PV = SV(NC,1> - OA 

IF (f)VI B50, 844, 850 

844 DT = DA/OA 
GO TO 960 

846 OV = SV( NC,K+l J - SV(NC,K) 

IF C OV » 85C, 848, 850 

843 DT = DA/SV(NC,K) 

GO TO 860 

B50 OT = (SVINC.K + ll - SORT ( SV{ NC ,K + 1 )* * 2 - DV+OA/C< 361) ) /(0V/DI14H 
860 CONTINUE 

TAUIII = C( 35) - $V 1 18 , K + l) - DT 
870 CONTINUE 
RF TURN 

890 FORMAT 15X2 7HSV ARRAY EXCEEDED ON CVC-E 12) 

END 

SUBROUTINE RC12(URZ4,URZ5) 

INTEGER UP Z4 

DIMENSION C (401 , XN(40,2), SV(L8,50), VTRANt 18) , RK(14,4), 

1 SC.U4,5), DC ( 1 4 1 , 4(6), 0(15), F(20), VOLUO), RMT<2), 

2 BC(4 } , 0F(6 I , TAU (5 ) , CC t 3 » , CH8(3), CH(4), CPH(3), 

3 00(41 

COHMON/Z/ C, XN, SV, VTRAN, RK, SC, Xi A, 0, F, VOL, RUT, 6C, QF, 

1 TAU, CC, CHB, CH, CPH, DQ, VE? VI* CPB, CPT, CADK, X, OT, 

2 IRK, LOC, I TE° X, INDEX, I, J ■ 8, N 

CU’«MON/R/ XOS.XHH.CXT , WORK, DUH l, 0UH2, 3UH3,W0RK2,RHTB, RHT82 ,TIHEOF 
I ,R«LIN 
CXT=C(35)+X0S-1C. 

IF(CXT.LE.0.)CXT=*0. 

C RESPIRATORY FPFOIJFNCY. 

FRF0=3.t + f.3I5*lPPTI2)TC(26)) 

C OE AD SPACE VENTILATION 

.OEAOVT=. 1 tOT^FREO* . OT85*VE 
C( jn = 10EADVT*C( 1)+VE*CU*U)/(D6ADVT+VE) 

C( 32) = (PC ADVT 4C ( 2 1+VC#0UM2 ) /( OEADVT +7E I 
C( 331= (DEAOVT*C( 3 ) +VF40UM3) /( DEADVT+VE 1 
MINUTE VCLU”F. 

TVNT = 0EADVTMVE+VI)/2. 

C HEART RATE. 

HR ATE=43. 8* (RMT (2 ) EC (26 1 1+54.5 
C 
C 

C SEE IF WORKLOAD HAS CHANGED. 


I 

iv 

I 



IFIURZ5 .NF. WORK) GO TO 500 
C 
C 

IFIMARKET.E0.0I GOTO 101 
l WURK=W0HK2 
RARKER = 1 
C 

C SYSTEM RESPONSES: TIME CONSTANTS FOR WORT LOAD LEVELSUNCREASING) . 
IF I WQPK.LE. 0. I GOT 02 
IFtwrm.GE.50. 1 T C T = 2 . 3/( 2.*Wf)RK/20C. 1 
IF IWORK.LF .*50 ) TCT = A *6 
C TISSUE 02 MBTA.IOUC RATE, 

C RMT12) COMFS FROM GUYTON. 

VTtME = l.l-I.l*EXPt-TCT#(CXT-T IMFCNI/1- 921 
C TERM USED IN VI THAT IS A COMPONENT OF It ANSIENT RESPONSE RELATED 
C TO WORK LOAD. 

RMLIN =SS02WIWPRX)-ISSOZ«(WORK)-RMTB2 *IW-VT!KE1 
IFIVT IHE.GE .1.1 RML 1N= S S()2W( WORK ) 

C TISSUF CC2 METAPTLIC RATE. 

PITH 80*RMT(2I 

IFITVNT.GT.37. I R MT <11 = I TVNT+A0.771 *FIIT( 21 /S8.5 
1F(URZA .EO. 0) GO TO 1230 
WRITE (6,333) RMTtl) ,RMT«2) 

333 FORMAT! 'O', IX, 25HCHANGE IN METABOLIC RATES, SX,7HMRC02=- ,F10.A, 
1 5X,6HMRQ2= ,F1C.A,/1 
2 CONTINUE 

c 

IF t UR ZA .£0. 01 GO TC 1230 
C 

C HERE IF GCING TO PP1NT0UT. 

C ARTFPIAL N2 tension. 

1210 PAN2 = 0(1)90(3) 

C TISSUF 02 TENSION. 

PT')2 =» C < 8) /Dt 3) 

C TISSUE N? TENSION. 

PIN2 = C (9) /0(A) 

C CEREBROSPINAL FUIIC PH , EQUATION 6.2 . 

PMSCF o q. - PCFIICH(A)) 

C VENOUS BRAIN H* CONCT NT RAT I ON , EQUATION A.T . 

HV6 » CADK*F(A)/(CC(2) - F < A* > 

C VENOUS BRAIN' PH , FOUAT ION A. 6 . 

PIIVO <■ 9. - PCFl(HVP) 

C VENOUS TISSUF F+ CONCENTRATION , EQUATION 5.7 . 

HVT * C.\nK*FI6) /(CCI3) - F ( 6 ) ) 

C VENOUS Ttssur PH , EQUATION 5.6 , 

PHVT » 9. - RCFltHVT) 

C RFSP! RATORY QUOTIENT (ALVEOLAR). 

RQ * 1 (C( L t )*VTP A M ( A ) + QF ( l ) +VTRAN 1 1 1 1 /C 1 10) - CCU))/ 
l (F(1) - (CUl)*VTRAN(5) * QF(1 J *VTRAN( B] l/CI 10) 1 

QF ( 5 ) * OF (6) - RO 

c 

WRITS (6,1510) CXT , RQ, QF{5) 

C 

WRITE (6,inl5) (Clli; I « 1,3), («(!), 1 « 1,31, F(7), FCl), 
l PA*I2 

WRITE (A, 1820) CC(l), F(9), FUO). F(7), Fill, PAW2 , CHCl), 
l CPH(l), CHB( l ) 

WRITE 1 6 , 18251 (CIO, I « A, 6), (C<,(1), I « A, 6), CPO, F(17), 

1 FUR), CHIP), CPH(2> 

WRITE (6 , L 8301 CCIM, I « 7,9), (DCIll, I « 7,91, CPT, PT02, 


I 

— Q 
U) 

I 



no n on o n o 


I PTN2,.CHI3), CPH (3) 


W'MTF 

(6,t8'35) 

tOCII), I = 

12,151, 

:cu), t = 

12, 1A) 

, CHIA), 

1 

PHCSF 






WRITE 

(6,1850) 

CC I 2 ) , FI 12 ) 

, C 1 6) , 

:pb, fut) 

, File) 

, Hve, 

I 

PH VM » CHfi (2 ) 





WRITE 

(6,1855) 

CC ( 3) , F I 13) 

, C<9), 

SPT, PT02, 

PTN2, 

HVT, 


l PHVT, CHfi (3 ) 

WRITE 16,18501 ITAUU), I = 1*5) , VI, VE, C(10), Ctll), 0CI101, 

1 DC (1 1 ) 

WRITE (0,1855) FREC,TVNT,OEADVT * HR ATE 

1230 RETURN ' 

1200 FU«»AT (SH XXX X 5X7F 1 0. A ) 

1202 FORMAT (RFIO.A) 

1805 FORMAT ( 1H1) 

1810 FORMAT I LHOOXAFTIMrFlO.A, TAX oHALV RO : 10 . A , 3X7HRQ 0IFF,F8.4/ 

1 16X3HCO’eX2HD2nx2HN27X21H) PRIVATIVE S9X4HPC026X 

2 3HPn2 7X3HPN2 7X4H<HM7X2HPH5X4HH)02> 

1815 FORMAT (3X8MALVF0LAR9F10.AI 

1870 FORMAT (3<“MARTFRIAL3F10.A,30X,5F10.A,Ffl.AI 
1825 FORMAT 1 6X5HRRA I NX l F 10 .4 ) 

1830 FORMAT ( 5X6HTI SSUE It F 1 0. A I 
1835 FORMAT 18X3HCSF3CX8F10.A) 

1 RAO FORMAT (AX7HV PP A IN3 F 10 . A, 30X , 5F 10. A . - 0 . A ) 

19A5 FORMAT (3X8HV T I SSUF3F10 . A ,30X , 5F10 . * ,F8. A ) 

1850 FORMAT I 5X1SHTR AN SPORT TIMES — AX2HA J 3X2H VB8X2HVT8X2HAT8X2HAC2X 
1 2H**AX2HV1SX2MVF3XIH09X2HFB7X11H0ER1VATIVES/21X, l OF 10.A»F8.A) 
1855 FORMA T { 3X , 9HRE SP FREC , F 10.4 , 5X , 13HM 1 N JT E VOLUME ,F10 .A ,5X, 

1 15HDEA0 SPACE VENT.F10.4, 10X,10HHfcAU RATE ,F10.A) 


Hc“E GUYTON SENT WORK LOAD CHANGE. 

500 »0RK2 = IJTZ5 

WRITE(6,305) H0RK2.CXT 

305 FORMAT 1*0*, 20 I * ^ * 1 / * CHG.1N WORK LOAD, HORKb 1 ,F7.2 , 'WATTS •, 
l MAT ',F10.A,*MINS)M 
TIMEDN^CXT 

SYSTEM RF.5P0NSFS: TI^E CONSTANTS FOR WOPS LOADS AND TISSUE 02 
METABOLIC RATF. 

ir (wnMK2.GE.W0R'< >RMTB2=RMTI2) 

DECREASING WORK LOADS 

IF I W0RK2.LT. WURK) RMIM=RHT<2) 

IF! RURR7.LT .WORK I RMT 8= SS02W ( WORK 2) 

IF UWDRK2.LT.WCRK) .ANU. [WORK. GE. 50. ) ) TCT»2 .3/ 12.*WDRK/200.1 
IF I HORK 2, LT. WORK), AND. (WORK. LT. 50.1) fCT*4.6 
I F f fORK2 . Gc .WORK ) G0T01 
101 W.iPKa^UP,K2 
MAP K FR* C 

TISSUE 02 MFTA'JC'LIC PATE. 

R U T(2) OXI'FS FRCM GUYTCN. 

VT !MF*l.l-l.mXP(-TCT*lCXT-T!KEnN)/3,34) 

C TEPM usfn IN V! THAT IS A COMPONENT CF TRANSIENT RESPONSE REIATEO 
C TO WORK. LOAD. 

RMLIN »RMTti-(RMTe-RMT'()«( 1,-VTJME) 

IFtVTlMF.OFd.) RMLIK.RNTa 
C TISSUE C02 METABOLIC PATE. 

. RMTUl = .88*RMTt21 

IF! TVNT.GT. 37. ) RMT ( 1 ) = ( TVNT+A0.77) *1 MT 12 I/8B.5 
IFUJRZA .EO. 0) GO TO 1230 
WRITE (6,333) RKT 1 1 ) ,RHT (21 
GO TO 2 


-■tlZ- 



o o 


END 

subroutine rc3 

OIMENSIUN C( 40 ), XN ( 40 , 2 ) , SV( 13 , 50 ), VTRANtlB), RKU 4 , 4 >, 

1 SCU4,5), CCI14), A I 6 ) * OU>D, F (20 ) » VCILUO), RMTI2), 

2 nCI4), QF(6>, TAU(5), CCC3). CHE(3), CH(4I, CPH(3>, 

3 GO 14 I 

CCIP.MON/Z/ C, XN, SV, VTRAN, RK, SC, DC, A, 0 , F, VOL, RMT , BC, OF, 

1 TAIJ, CC, CHB, CH, CPH, 00 , VE, VI, CPB, CPT, CAOK, X, OT, 

2 IRK, LCC , 1 TERX, INDEX , I, J, . 1 , N 

C 6969 FORMATUH 7 HSUB RC 3 ) 

SETS T I PF-DEP ENGENT EXPRESSIONS 
TISSUE BLOOD ELPO. 

of i i i = cuoi - cun 

C ARTEUAI C 2 TRASKS. 

Ft 1) = 0111*012) 

C ARTERIAL CG 2 CONCENT RAT ION , 

fi?) = Dtsncm 

C BRAIN 02 CIjNC P i NTR AI ON / ( CONV.FAC TnR* 9 Gl')BI LI TY CCEFF.FUR 02 ) 

F t 3 ) = C( 5 >/D( 3 ) 

C IC 0 NV.FACT()P*$ 0 LIJ 6 ILITV COEFF.FOR CD 2 ) >■ BRAIN C 02 TENSION. 

F( 4 > = D ( 2 ) *CP fl 

C TISSUE 72 CGNCFN'T R AT I CN / (CCNV. FACTCRA 33 LUB IL ITY COEFF.FOP 02 ) 

FI 5 ) = C ( 3 ) /Dt 3 ) 

C ICCNV.FACTCR*SCLUEILITY COEFF.FOR C 02 ) * TISSUE C 02 TENSION. 

F I 6 ) = ni 2 )*CPT 
C ARTERIAL CO? TENSION. 

FIT) = 0 |I)*CI 1 'I 
C ARTERIAL 02 TENSION. 

FIRI = Oil ) *C( 2 ) 

RETURN 

ENO 

SUBROUTINE RCA 

DIMENSION C( 40 l , XN( 40 , 2 ), SVU 8 .S 0 ). V)RAN(IB), RK 1 14 , 4 ) , 

1 SCI I A, 5 ) , DCI 1 AI, A 1 6 ) , Dili), F( 20 ), VOLIIO), RMTI 2 ) , 

2 BC(A), OF 1 6 I , TAU ( 5 ) , CCI 3 '. CHB 13 ), CH(A), CPHI 3 ), 

3 00 (A) 

COMMON/?/ C, XN, S V, VTRAN, RK, SC, Cl, A, D, F, VOL, RMT, BC, OF, 

1 TAN, CC, CMB , CH, CPH, DO, VE. VI, CPB, CPT," CAOK, X, DT, 

2 IRK, LUC, ITERX, INDEX, I, J, K, N 

C ITERATES FOR CCtll, ARTFRIAL C 02 CONCENTRATION 
C 6969 FORMAT ( IH 7 HSUB RCA) 

AlO CALL RC 2 1 (CHB(L), Fill, F( 2 ), CCtl), CHU), CPH(D) 

X = ICCIl) - F( 2 )J/( 0 . 01 *P( 7 ) ) 

X = RCFt(X) 

C SEE EOUATICN 3 . 1 , X= CA 1 C 02 ) . 

X = HC(l) + 0 . 375 * (CI 17 ) - CHB ( 1 ) ) ♦ = 12 ) - D 18 )*(X - 0 . 14 ) 

C CCll) * CA ( C 02 ) . 

CALL RC6 (CC ( L ) ) 

CCU) = CCtl) ♦ 2 . 0 * ( X - CC ( l ) ) / 3.0 
C 3000 PURMATUH , 5 HCC( 11 , 5 X,E 16 . 6 ) 

' IF (. 1 TERX) A 20 , AlO, A 20 
A 20 RETURN 
END 

SUBROUTINE RC 5 (CP, Ffl, CCB, BHC ) 

DIMENSION C( AO), XN( 40 , 2 ) , SV( 18 , 50 >, VTRANIIB), RKdA.A), 

1 SC( 14 , 5 ) , CCUA), A ( 6 ) , 0 ()>), F( 20 ), VOLIIO), RMT (2 ) , 

2 6 C( A ) , 0 FI 6 ), TAUI 5 ), CCI 3 ), CHB( 3 ), CH(A), CPHI 3 ) , 


-0 

U) 

I 



oonn 


3 no(A) 

CHMMON/Z/ c, X N, SV, VTRAN, RK, SC. lit., A, 0, F, VCL, RMT, BC, OF. 

1 TAU, CC. CHH, CM. CRH, DP. VE, VI. CPB, CRT, CAOK, X. DT, 

2 IRK, Lac, ITPRX, INOEX, I, J, II, N 

C IT SPATES FPR BRAIN AND TISSUE PC02 
C6569 FOHPATUH 7HSUB RC5I 

510 X = ( CCD - FB1/(0.01*CP! 

X = RCFItX) 

C SEE EOUATION 4.1, X = PBICr2) . 

X = (-BHC ♦, CCB * 0(i0)«(X - 0.14D/U 2) 

C CP' PBICC2 I . 

CALL RC6 ICPI 

CP = CP ♦ (X - CPI/10.0 

c cerebral »Lcrn flcw. 

FB = 0(2) *CP 

C3000 F0° MAT (1H ,4HCP = ,E16.6,4HFB= E16.6,SHCCB= El 6. 6, 5HftHC = E16.6) 

IF UTERX) 520. 510, 520 
52 0 RETURN 
END 

SUBROUTINE RCft t Y> 

D1PFNSICN C ( AO ) , XN ( 40 , 2 1 » SVU3.50), VTR ANt 1 8) , RK ( 1 A, 4) , 

1 SC( 14,5) , DC ( 1 4 I , A(6), Dt L!, , F120), VOLUO), RMTI2I, 

2 BC(4), 0 F [ 6 ) , TAU15I, CCt 31 CHBI 3) , CH(4), CPril 3 1 , 

3 00(4) 

COWON /It C, XN, SV, VTRAN, RX , SC, U, A, D, F, VOL, RMT, BC, OF, 

1 TAU, CC, CHB, CH, CPH, tip, VE, VI, CPB, CPT, CAOK, X, DT, 

2 IRK, LOC, ITFRX, INDEX, I, J, H, N 
CHECKS CONVERGENCE OF ITERATIVE PROCEDURES 

RC4 : X-CA l C02 I , Y*CCU1 . 

RC5 S X»PB(C02), V-CP . 

KC19 : X=CVB(C(I2>, V=CVC . 

C6969 F0R m AT(1H THSUB RC6) 

ITERX = 0 

IMFF = A6S (IX - YI/Y'J 
IF (OIFF - I.OF-5) 620, 620, 630 

620 ITERX = 1 
630 RETURN 
END 

SUBROUTINE KC7 

DIMENSION C(40), XN(40,2I, SVII8.50), VTRAN1 18 1 , RK(14,4), 

1 SCt 14,5 I , DCt 14 I , At 6 1 , D( 15) , F120), VOLUO), RMT12) , 

2 BC 1 4 ) , OF ( 6 ) , T AUt 5 ) , CC 1 3 ) , (HB[3), CH(4), CPH(3l, 

3 00(4) 

COWON/ZA C, XN, SV, VTRAN, RK, SC, t(.. A, D, F, VOL, RMT, BC, OF, 

1 TAU, CC, CHB, CH, CPH, DO, VE, VI, CPB, CPT, CAOK, X, OT, 

2 I P K, LOC, ITERX, INDEX, I, J, II, N 

Cn“MONAPA XOSiXMHfCXT, WORK,DUH l ,0UM2 tUUM3, W0RK2.RMT8, RMTB2.TJME0F 
1 , R M L IN 

C6969 FORMAT ( 1H 7HSUB RC7I 

C FILLS SV ARRAY WITH INITIAL CONDITION 

CALL RC16 

IFtXDS.GT.XHHI GPT02 
DO 725 I * 1,17 

DO 720 J = 2,50 

SVU,J) = SVU.l) 

720 CONTINUE . 

725 CONTINUE 

2 CONTINUE 


- 76 - 



on 730 J a 2,50 

SVI 13 , J) = SV(1B,J - n - 0(141 
730 CONTINUE 

C3000 F0P4ATUH , 12H18SV S 0 ( 14 1 , 6{ 3X, E16. 6) t IH ' ,6 ( 3X, E16.6) /IK ,7I3X,E1 
C C6.611 

RETURN 
ENC 

SUBROUTINE PC9 

DIMENSION C( 40 1 « XNt40,2), SVU8,50), VTRAN( 18) , RK(14,4), 

1 'SCU4,5J, cr,(l4). A (6 ) t DI15I. F ( 201 , VOLdOl, RMTI2I, 

2 BC 1 4 ) » OF ( 6) r TAU( 5) , CC(3), CHH(3), CH{4), CPH(31» 

3 00(41 

COWPON/Z/ C, XN, SV, VTRAN, RK, SC, CC, A, 0, F, VOL, RHT , 8C, OF, 

1 TAU, CC, CUB , CH, CPH, 00, VE , VI, CPB, CPT, CAOK, X, DT, 

2 IRK, IOC, ITSRX, INDEX, I, J, M, N 

C SETS VALUES IN VTRAN ARRAY 

C6969 FORMATUH 7HSU0 RC 9 ) 

on 960 l = 1,5 

TA = T AU ( 1 1 ~ (C(35) - SVf 18,1)1 
LOC = TA/ 0 ( 14 1 

IF ( LEC - 49J 904, 904, 902 

902 WRITE (6,990) I, LOC 
IOC = 49 
5C4 XLOC = LOC 

T 8 = XLOCRO ( 14 1 
OT = TA - TB 

GO TO (910,920,930,940,950), ( 

910 00 914 J = 1,3 

C LUNG TO BRAIN C 02 ,r 2 ,N 2 TIME DELAYED ARTERIAL CONCENTRATIONS. 

VTRA’-I (J) = RCF3 ( J ) 

914 CONTINUE 

C LONG TO BRAIN Hf TIME DELAYED ARTERIAL CONCENTRATION. 

VTPANU5) » RCF3{13) 

GO TO 960 

920 DO 924 J = 4,6' 

C BRAIN TO LUNG C02,02,N2 TIME DELAYEO VEIf 3US CONCENTRATIONS. 

VTRAN I J 1 = RCF31J) 

924 CONTINUE 

C BRAIN TO LUNG COMBINED C02,02,N2 TIME DE..AYED VFNOUS CONCENTRATIONS. 
VTPAN ( 16) = RCF3I16) 

GO TO 960 

939 00 934 J = 7,9 

C TI5S1JF T) (UNG C02.C2.N2 TIME DELAYED VENOUS CONCENTRATIONS. 
vtpanij) = or.F3(JI 
934 CONTINUE 

C TISSUE TO LUNG COMBINED C02,02,N2 TIME DELAYED VENOUS CONCENTRATIONS. 
VTRANU7) = PCF3U71 
GO TO 960 

940 CO 944 J = 1,3 

C LUNG TO TISSUE C02 ,02 ,N2 TIME DELAYEO ARTERIAL CONCENTRATIONS. 

VTRAN I J*9 I = RCF3U) 

944 CCNTiNU? 

G 1 TO 960 

C LUNG TO CAROTID SITE H* TIME DELAYED ARTERIAL CONCENTRATION. 

950 VTPAN ( 13 ) = PCF3I131 

C LUNG TO CAROTID SITE 02 TIME DELAYED ARTERIAL TENSION. 

VTRAN ( 14 1 * RCF3I 14) 

960 CONTINUE 

C NAMEL I ST /OONM/VTRAN 


1 

-J 

I 



RETURN 

990 FORMAT ( 5X27HSV ARRAY EXCEEDED CN CYCLE 1 2 » 12H WITH LOC = ,14) 

ENO 

SUBROUTINE RCIO 

0 1 M ENS ICN C (40 ) , XN( 40 , 2 ) , SV(L8,50), VTRANU8), RKU4,4), 

1 SC(14,5), CC ( 14) , 4(6), 0(13), F(20), VUL(IO), RMT { 2 1 ? 

2 BC 1 4 } , 0F( 6 ) , TAU(S), CC(3), CH8(3), CH(4>, CPH(3), 

3 00(4) 

COHMON/Z/ C, XN, SV, VT RAN , RK , SC, CC, A, I), F, VCL, RMT, BC, OF, 

1 TAU, CC, CHH, CH, CPH, DQ, Vfc, VI, CPB, CPT, CAOK, X, DT, 

2 ( R K , LOC, ITCHX, INDEX, I, J, R, N 

C6969 FORMAT! 1H BHSUB RCIO) 

C COMPUTES EMPIRICAL FUNCTIONS FOR ACR3IAC OUTPUT AND BRAIN BLOOD 
C FLOW DIFFERENTIAL EQUATIONS 

00(1) = C. 

00(2) = 0. 

C Fiei = PA ( 02) . 

IF (F(S) - 104.0) 100B, 1G24, L024 

C ( DELTA ) CB ( C2 ) , EQUATION 7.9 . 

100B 00(2) = II < 7.65S9E-8AF (8) - 2. 324E-5>AF (8) ♦ 2.6032E-3 ) “FIB ) 
t - 0. 1323 )*F (a ) 4 2.785 

1016 )F (00(2)) 1024, 1044, 1044 

1024 CQ (2 ) = 0.0 
C 

C IF PCC2 GT 60 OQ (3 1 STAYS AT ITS VAL^E AT 60 OLO ROUTINE SETS 

C THE VALUE OF 00(3) EQUAL TO 0 
C 

1044 IF IF ( 7) - 38.0) 1043, 1052, 1052 

c toELTAioetcnai , equation ?.u . 

1043 00(4) = (B.0163E— 4FFI7) - 3. 1073E-2 ) *F ( 7) 4 2.3232E-2 
RETURN 

1052 IF ( F (7 ) - 44.0) 1056, 1016, 1060 

1056 00(4) * 0.0 
RETURN 

C t CELT A ) COt CC2 ) , EQUATION 7.13 . 

1060 r,n ( 4 » =■ ( ( (—2. 1 74 BE- 7* F ( 7) 4 9.391BE-5)*F(7) - 1.2947E-2) *F { 7) 

1 4 C.76C7)*F(7) - 15,58 

C AAPELIST/C6/CQ.F 

RFTURN 
END 

SUBROUTINE RCI l 

01 MENS ION C (40 ) , XN(40,2), SV(1IS,50). VTRANS18), RK(14,4), 

1 SCI 14,5), DC 1 1 4 I « A (6 ) , 0(15) , F120), VOL(IO), RMTI2), 

2 BC ( 4 ) , <3F(6), TAU(5), CCI3I, CHBI3), CH(4), CPH(3), 

3 DQ(4) 

COMMON/Z/ C, XN, SV, VTPAN, RK, SC, (C, A, D, F, VOL, RMT, BC, QF, 

1 TAU, CC, CUB, CH, CPH, DO, VE . VI, CPB, CPT, CADK, X, DT, 

2 IRK, LOC, ! TERX , INDEX, I, J, R, N 

C CALCULATES DIFFERENTIAL FQUATIONS 

C6969 FORMAT! LH BHSUB HCll) 

CALL* RCI 7 
C EQUATION 10. I . 

Dell) *> (VI«C(3l) - VE*C(l) 4 D(9)*(f 1 11)*VTRAN(4) 4 QFIl) 

1 RVTRAflt 71 - C(10I*CC(I1))/C{>2) 

C EQUATION 10.2 . 

• DC 12) ■ ( Vl*C( 32) - VE*C(2) 4 D(9 >*<i 1) 1)*VTRAN( 51 t QFU) 

I OVTPANIS) - C(lO)*F(9m/Ct2. r :) 

C EQUATION 10.3 . 


I 

-3 

03 

I 


0CI3) = ( VI *Ct 33) - VE*C (3) + 0 ( 9 I * ( C I 1 1 ) * VTRAN C6) «■ QFU1 
1 *VTRANI9> - C(l0J*F{10)))/C(2iJ 

C equation 10.4 . 

OC 1 4 ) = ICI25) ♦ CtU)MVTRAN(l) - CC12II - F(14))/C(23) 

C EQUATION 10. b . 

DC ( 5) = 1-CI26I + C(11)*IVTRAN12) - FI12I) - F(15))/C(23) 

C EOUATION 10.6 . 

DClbJ s I C( 11 ) ♦(VTRAN (3) - C(6)J - F ( I 6 ) ) /C ( 2 3 ) 

C EQUATION 10.7 . 

DC (7 > = (RMTU) ♦ OF 1 1 >*( VTRAN ( 10) - C C ( 3) I ) /Cl 24) 

C EQUATION 10.8 . 

OC ( 3 ) = 1-PHTC2) R QF ( l )* I VTR AN( 1 1 ) - F 1 13 ) ) ) /C 1 24 ) 

C EOUATICN 10.9 . 

OC l 9) - 0Fll)*(VTPAN(12) - C(9))/C(24) 

OC 1 10 ) = 0. 

C 

C 

C OEPENDANCE OF CARDIAC OUTPUT ON TISSUE 
C UTILIZATION OF OXYGEN. 

XA8- 5. 5 * (RMTI2I-.21 5I+6.-CI 10) 

tFt(RMT(2).GT..2L5).AN0.(XAB.GT.Q.))0t(l0)*DCU01«-XAB/.0l0 

C 

C 

C EQUATION 7.7 . 

octlt) = l-Clll) + 0.75 ♦ 0012) ♦ 0QH))/CU9) 

C FQUATIUN 1.10 . 

net 12 ) = F( 141/101341*0(11)) 

C EQUATION l.ll . 

DC ( 1 3 1 ? Ft 151/IC134)*0( 12)) 

C EQUATION’ 1.12 . 

nc(14) = F(.16)/(C(34I»C<13I) 

C NA«FL IST/AB/DC 
RETURN 
FNO 

SUBROUTINE RC13 

01 v FNSION C. ( 4 0 ) , XN(40,2), SV(18,50), VTRAN11B), RKI14,4), 

1 50(14, 5), DC 1 141 , A ( 6 1 , 0051, F (20) , VOCttO), RHT (2 ) , 

2 BCI4), OF 16 1, TAU(5), CCI3), CHQ(3), CH(4), CPHl 3) » 

3 00(4) 

COMMON/2/ c, XN, SV, VTRAN, RK, SC, fit. , A, 0, F, VOL, RMT , BC, QF, 

1 TAU, CC, CHD, CHi CPH, PQ, VE.-VI, CPB, CPT, CAOK, X, OT, 

2 IRK, me, ITERX, INDEX, I, J, II, N 

C6969 . FORMAT (1H BHSUB RC13) 

C SOLVES M DIFFERENTIAL EQUATIONS BY FOURTH-ORDER RUNGE-KUTTA AND 

C ACAMS-MPUTLFN PREDICTOR-CORRECTOR METHODS 

C tvAMELlST/OP(j/C,CC,SC 

’ IF IISK - 41 1304, 1356, 1356 

1304 03 1362 INDEX = 1,4 

, PD 1309 I ■= l,M 

RK( I, INDEX) * DCIl) 

1303 CONTINUE 

GO TO (1312, 1320, 1328, 1340), INDEX 
1312 CO 1316 I = l,M 

SC( I , IRK+1 ) » cm 
SCll, IRK) a DC 1 1 1 
1316 CONTINUE 
TI = CI35) 

13P0 COS) “ T! *• CO6I/2.0 
DO 1324 1 * l»* 


l 



CCIJ « SCtl,!RK+l) + C(36)*RKtI, INDEX *2.0 
1324 CONTINUE 
GH TO 1336 

1328 Ct 351 = T[ ♦ C [ 36) 

CO 1332 I = 1,N 

Ctl) * SC 1 1 * I RK*1 I + Ct36)*RKtl, INDEX* 

1332 CONTINUE 
1336 CALL RC14 
. GO TO 1352 
1340 00 1344 I = it M 

CU) = SCH,tRK+ll + C(36)*tRKU,l) i 2.0*RKtI,2) * 2.0*RKtI,3> 

1 + RK (I , 4) I /6. C 

1344 CONTINUE 

IRK = IRK ♦ 1 
1352 CONTINUE 
RETURN 

1356 DO 1360 l » 1,M 
SCII, 5) = Ct I) 

sc ( 1 , 41 = net n 

cm » SC 1 1 r 5 1 + C t 36) * ( 55. 0*SC 1 1 ,4) " 59.0*SC(I,3) «■ 37.0«$C(I,21 
1 - 9.09SCI 1,1) 1/24.0 

1360 CONTINUE 

Ct 35) = Ct 35 1 * C t 36) 

NC35=Ct35)/C( 361 + .1 
C ( 151 =C ( 36) *NC35 
1364 CALL RC 14 

OO 1363 I = 1,M 
SC (1 , 1) = C 1 1 1 

CII) = sen, 5) ♦ Ct36)*t9.0*DCtn ♦ 19,0* SCI 1 ,4) - 5.0* SCI 1 ,3) 

1 * sell, 211/24.0 

1363 CONTINUE 

DO 1372 I = 1,M 

IF (A3S ten) - set 1 , II ) - l.OE-3) *372, 1372, 1364 

1372 CONTINUE 

no 1376 I « 1,H 
DO 1376 J = 1,3 
SC(I,J) = SC(I,J+l) 

1376 CONTINUE 
RFTURN 
END 

SUBROUTINE RCt 4 

DIMENSION CI40), XNI40.2), SV(18,50), VTRANt 18), RKU4,4I, 

1 SCt 14,5 ) • DC 1 14 ) , A<6), 0(16), F(20), VOLtlO), RHTI2), 

2 DC ( 4 ) , OF 1 6 ) , T4UIS), CCI31- CHBI3), CHt4), C PH 1 3 ) , 

3 PDC4) 

C01R0N/Z/ C, XN, SV, VTKAN, RK, SC, tL, A, 0, F, VOL, RMT, BC, QF, 

1 TAU, CC, CHB, CH, CRH, DO, VE, VI, CPB, CPT, CADK, X, DT, 

2 IRK, LtJC, ITERX, INDEX, I, J, H, N 

C CALLS OTHER SUBROUTINES IN C BLOCK 
C&969 FORMAT! 1H 8HSUB RC14) 

CALL RC3 
CALL ‘RC8 
CALL RC9 
CALL RC4 

CALL RC 5 (CP8, F(4), Ct 4 1 , BCt2>) 

CALL RC21 ( CHB 12), F<3), F(4>, C(4), (Ht2), CPHt2)l 
CALL RC 19 (CPB, CHB ( 2 ) , CCt2), BCtU, Ft4)) 

CALL RC5 tCPT , F{6I, C( 7), t)C(3JI 

CALL RC21 1 CHB 13) » ftS), Ft6), Ct7), I.HOI* CPHOII 



o o 


CALL RC 1 9 (CPT, CHB(3), CC(3), BC(11, FI6I1 

CALL RC 10 

CALL RC20 

CALL RCI 1 

RETURN 

ENO 

SUBROUTINE PC15 

DIHFNSION C(4QJ, XNI40.2), SV(18,50), VTRAM161, RK114.4J, 

1 SCI 14,51, CCI14), A ( 6 I » D(15), FI20J, VOLllOJ, RMTI2) , 

2 BC(4), QF(6>, T AUI5 1 1 CC(3>. CH0I3I, CHI4*, CPHI31, 

• 3 00(4) 

COPMON/Z/ C, XN, SV, VTRAN, RK, SC, CK , A, D, F, VOL, RMT, BC» OF, 
1 TAU, CC, CHS , CH, CPH, 00, VE, W, CPB, CPT, CAOK, X, OT, 

Z IRK, LOC, 5TERX, INDEX, I, J, A , N 

CA969 FORMAT ( 1H 8HSUB RCI5J 

C NAMELIST/SCH/SV 

C SHIFTS VALUES IN SV ARRAV 

00 1530 I = 1,18 
00 1520 J = X, 49 
JM = 51 - J 
JMM = JM - 1 
SV{ 1 , JK ) = SV( I ,JMH) 

1520 CONTINUE 
1530 CONTINUE 
RETURN 
END 

SUBROUTINE RC16 

DIMENSION C 140 1 , XNI40 , 2 1 1 SV(X8,50I, VTR AN (181, RK(X4,4I, 

1 SCI 14, 51 , DC (141, AI6) , D(15), FI20), VOLIlOl, RHTI2I, 

2 PCI4) , OF 1 6 ) , TAUI5), CCI3), CHBI3), CH(4), CPH13), 

3 00(4) 

COMMON /Z/ C, XN, SV, VTRAN, RK, SC, DC, A, D, F, VOL, RMT, BC» OF, 

1 TAU, CC, CHB, CH, CPH, OO, VE, VI, CPB, CPT, CADK, X, DT, 

2 IRK, LOC, ITFRX, INDEX, I, J, t, N 

COMMON/ P/ XDS,X«H,CXT, WORK.nUMl, D0M2,f UM3,W0RK2,RMT8,RNTB2,TIH£0F 
X *RPL I N 

C6969 FORMAT! IH 8HSUB RC16) 

SETS VALUES FOR SV ARRAY 
ARTFRIAL C02 CCMC ENTPAT ION. 

sv<t,n = cci u 

C ARTERIAL C2 CONCENTRATION. 

SV(2,1) = F(9) 

C BRAIN VE ACUS Cr? COAXENTRAT ION. 

SV(4, II = CC<2) 

C ARTFRIAL N2 CONCENTRATION. 

SVI3.1) = F I I 0 ) 

C BRAIN VENOUS 0? CONCENTRATION 
SVI5, 11 = FI121 

C PR AIN VEACUS A2 CCNCENTRAT ICN. 

SV ( 6 , 1 1 = C(6) 

C TISSUE VENOUS C02 CONCENTRATION. 

SV { 7 , l ) = CCI 31 

C TISSUE VENOUS 02 CONCENTRATION. 

SVI8.1J = F l 1 3 1 

C TISSUE VENOUS N2 CONCENTRATION. 

SVIR.U = C(9l 
C CARDIAC OUTPUT. 

, sviio.u « ciio) 

» 


I 

CO 

H 
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C CFRERRAL Bl CCD FLOW. 

svm.ti = can 

C TlSSUC BLCCD FLOW. 

SV112.D «= CF C 1 ) 

C flBTEMM. H+ CONCENTRATION. 

SVU3,1> = CMtll 
C ARTEP I OL OZ TENSION. 

SVH4.1) = FUJ 
C INITIAL TI“g. 

SVU5,1> = C.O 

C TOTAL GAS CONCENTRATIONS AT BRAIN EXIT. 

SVll6.il = SVl4.ll + SV(5,1) ♦ SVt6.il 
C TOTAL GAS CONCENTRATIONS AT TISSUE EXIT. 

sv c i t , 1 1 = sv(7,n ♦ sv(a.i) ♦ svtq.ii 

C SIMULATED TtME. 

SVtlB.l) = CI351 

RETURN 

END 

SUBROUTINE RC17 

0 1 MENS ION C (AO) i XN<40,2), SVIL8.S0I , VTRANI18I, RM14,4), 

1 SCI 14 1 5 > I DC t 14). AI6i, D(15>, FI20I, VOUIO), RMTI2I, 

2 RCI4). OF 16 ) . TAUI5J, CCI3>. CH8I3), CHI4I, CPHI3), 

3 0014) 

C0MP0N7Z/ C, XN, SV, VTRAN, RK, SC, CC, A, 0, F, VOL, RMT, BC, QF, 

1 TAU, CC, CH8 , CH, CRH, 00, VE, VI, CRB, CPT, CAOK, X, OT, 

2 IRK, LOC, ITEPX, INDEX, I, J, M* N 

COMMDN/R/ XOS, XMK.CXT ,NG8K, OU Ml » CIJM2 *DUM3 i WORK Z, RMT B, RMTB2,T IMEDF 
I ,RML IN 

C NAMELIST/BA C/ CHI4), CA OK ,C(ll),Ct 12 ),ECI4), Cl 37), 0(38) » VTRANI 14) » 

C 1TER«,VI,C(20) »C(16),VT°ANt15),C(2ll,\TRANtl3),C(37>,D(9),Ctll)» 

C ZVTHANtl6),0Ft l ), VTRANI l7),C(lO),F(lil, 

C6960 FOP»ATUH 8HSUB RC17) 

CALCULATES VENT 1 LATICN 
CFS H* CIINCC'ITKATION , EQUATION 6.1 . 

CH( 4 ) = C Af!K*C( ll)*Ctl2)/BC(4) 

IF ICI37) .GT. l.OE-5) GO TO 1768 

1704 VI = CI38) 

GO TO 1730 

1709 TERM = O.C 

OPC I S ION CN ARTFRIAL C2 TENSION AT CARL". ID BODIES'SITE. 

IF (VTRANtIM - 104.0) 1710, 1720, 1 T20 

1710 TERM = l2).6E-9]*t(104.0 - VT R AM 1 1 4 ) **4.9) 

CONTROLLER EQUATION AS A FUNCTION OF HI MORAL TERMS. 

1720 VI = CI20I*IC11AI*VTRAN< 15) * II .0 *• C t 16 ) ) *CH (4) > 

1 +, C (2L I *VTR An! 13 ) + TERM - C< 3"') 

INCLUSION OF NEURAL CC“PnN6NT AS A FUNC'ION OF HORK LOAD. 

5VNT2=SSVEN I ( SSQ2W ( WORK ) ) -VI 

irtlSVNT2.Gr.O.).AN0.tSVNT2.L6.l5.); VI»VI+SVNT2 
IFtSVNT2.GT.L5.) VI =VI + 15* 

DESCRIPTION CF 'transient VENTILATION RESPONSE. 

- SVNT = SSVENT I RML I N ) -VI 
I F ( SVNT.GT .0 .5 ) VI«VI*a.7S*SVNT 

EXPIRED VENTILATION RATE, EQUATION U.l . 

1730 VE = VI + DI9)*(CtlI)*VTRAN(l6) + QPU)*VTRANtl7) - CUO>*FtlU) 

IF (VI .LT, 0.0 -OR. VE .IT. 0.0) GO TO 1740 
RETURN 

1740 Vt » 0.0 



VE = 0.0 
RF TURN 
FND 


SUBROUTINE PC19 (CPA, CVHBA , CVC, BHCA, FC I 

DIMENSION C (AO) , XN(40,2), SV(IB,50>, VTRAN(IB), RK<14,4), 

1 SCI 14,5) , DC ( 14), A 1 6) f 0(15), F120), VCILUQ), RMT (2 ) » 

2 BCf4 ) , OF (6 ) i TAUI5), CC( 3 1 r CHB(3>, CH(4), CPH13), 

3 00(4) 

COMMON/ 2/ C, XN, SV, VTRAN, RK,' SC, DC, A, D, F, VOL, RMT, BC", OF, 

1 TAU>, CC, CHB, CH, CPH, 00, VE, VI, CPB, CPT, CADK , X, OT, 

2 IP.K, LUC, ITERX, INDEX, I , J, 1, N 
C NAHEL 1ST /DH2/CPA, CVHBA, CVC, BHCA,FC 

C 6969 FORMAT ( 1H 8HSUB RC19) 

C ITERATES FOR VENOUS eiAIN AND VENOUS TISSUE CD2 CONCENTRATION 
C TERM USED IN' EQUATION 4.2 . 

1910 X = (CVC - FC ) /( 0, 01 *CPA) 

C LOGARITHM SUBROUTINE. 

X = RCFl(X) 

C EQUATION 4.2 . 

X = RHCA ♦ 0.375* (C ( IT I - CVHBA) - 0!B)*U - 0.141 ♦ FC 
CALL RC6 (CVC) 

CVC = CVC + 2.0*(X - CVO/3.0 
IF (ITERX) 1920, 1910, 1920 

1920 CONTINUE 
RETURN 
END 


SUBROUTINE RC20 

01 MFNSI ON C (40) , XN ( 40 ,2 ) , SV(18,50), VTRAN(IB), RK114.4), 

1 SC( 14, 5) , DC ( 14 ) , A ( 6 ) , 0(151, F(20), VOL(IO), RHT (2 ) , 

2 BC (4 ) , OF ( 6 ), TAU( 5 ), CC(3), CHB(3), CH(4), CPHt3), 

3 no ( 4.) 

COMMON/Z/ C, XN, SV, VTRAN, RK, SC, DC, A, 0, F, VOL, RHT, BC, OF, 

1 TAU, CC, CHR, CH, CPH, OQ, VE, VI, CPB, CPT, CAW, X, OT, 

2 IRK, LOC, ITERX, INDEX, I, J, M, N 

C NAMELIST/NMF/F 

C6969 FORMAT (1 H BFSUB RC20I 

C SETS t IMF OtP EN CENT EXPRESSIONS 

C ARTERIAL OXYGEN CCNCENTRAT IUN INCLUDING EFFECTS OF HEMOGLOBIN. 

F(9) = DI6)«C(2) + CHStl) 

C ARTERIAL NITROGEN CONCENTRATION. 

Flld = n 17 ) *C » 31 

C TOTAL ARTERIAL GAS CONCENTRATION AT LUNG. EXIT. 

F(ll) = CC(1) + F( 9 ) + F(10) 

C VENOUS 9RAIN PXYGFN CONCENTRATION INCLUDING EFFECTS OF HEMOGLOBIN. 

FI 12) = C 15 1 ♦ CHB ( 2 ) 

C VENOUS ri'SSUC OXYGEN CONCENTRATION INCLUDING EFFECTS OF HEMOGLOBIN 
Ft 13) = C(8) + CHB ( 3 1 
C CXYGFN TENSION IN BRAIN. 

FI17) = C ( 5 ) /0 1 3 1 
C NITROGEN TENSION IN BRAIN. 

FU8I - C(6)/0(4) 

C PROCUCT OF DIFFUSION COEFFS.AND DIFFERENTIAL BRAIN - CSF GAS TENSIONS 
F ( 14) =» C(27)*t CP8 - C( 12) ) 

Ft 15) *> C (2 8 )R t F ( 17 > - C(13)) 

F( 16) » C129)*(FI 191 - C( 14) ) 


RETURN 

ENO 
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SUBROUTING 5 C2 1 (CHBA, Fit FC, CCA, (.HA, CPHA) 

DIMENSION C 140), XN 1 4 0 » 2 ) » SVII8.50) VTRANU8), RK [ 1 A, A) , 

1 SC ( LA, 5 I , CCI14), A 1 6 ) , 0(.5», FI20I, VOLIIO), RMTI2I, 

2 BCtAI, QF (6) , TAU 1 5 ) , CCI3 , CHB13), CH(A>, CPHI3), 

3 DO I A ) 

COM«ON/ZA C, XN, SV, VTRAN, RK, SC, IlC, A, 0, F, VOL, RMT, BC, QF, 

1 TAU, CC, CHR, CH, CPH, OQ, VI VI, CPB, CPT , CADK, X, DT, 

2 IRK, LOC, IT FAX , INDEX, I, J, H, N 

C6969 FCR^ATUH 8HSIIB RC21) 

C NAMEUST//>B/CHBA,FA,FD,CCA,CHA,CPHA 

C COMPUTES H+ ION, PH, AND OXYHEMQGLOt' N 

C ARTERIAL H, CONCENTRATION. 

CHA = CADK*FD/tCCA - FD) 

C ARTFPIAL PH. 

CPHA = 9.0 - RCFKCHA) 

C DEVELOPMENT OF EXPRESSION USED IN CALCULATION OF ARTERIAL 
C OXYHF W 0G L CM I N SATURATION. 

X = R CF2 t CPHA 1 
X = -X A FA 

x = i i.o - exp ixm *2 

X=ABS(X1 

c 

C AR1ERIAL OXYHEMOGLOBIN CONCENTRATION. 

CHfiA * X*Ctl7l 

RETURN 

END 

FUNCTION PCF1IH1 

DIMENSION CIAOI, XN(A0,2I, SVC IB, 501, VTRANI18), RKC1A.A), 

1 SCI I A ,5) , DC { 1 A I , A (6 ) , 01:5), FI20!, VDL(IO), RMTI2), 

2 PC I A ) , QF 161 , TAU 1 5) , CC131, CHBI3), CHtAl, CPHI3I, 

3 DO(A) 

COMMDN/Z/ C, XN, SV, VTRAN, RK, SC, DC, A, D, F, VOL, RMT, BC, QF, 

1 TAU, CC, CHE), CH, CPH, DQ, VE, VI, CPB, CPT, CAOK, X, DT, 

2 IRK, MIC, ITEPX, INDEX, I, J, M , N 

C LOGARITHM TO BASE 10 

RCF1 a 0.A3M29AA8 * ALOGIHI 

RFTURN 

END 

FUNCTION RCF2IZZ) 

DIMENSION CIAOI, XNIA0.21, SVC18,50‘. VTRAN(lfl) , RK11A.A), 

1 SCI 1 A, 5 1 , DC 1 1 A), A 1 6 1 , D> 15) , F(20], VOLIIO), RMTI2I, 

2 BCIA), QF<6), TAU(S), CC C il* CH8I3), CHCAl, CPHI31, 

3 DO I A I 

COMFON/Z/ C, XN, SV, VTRAN, RK, SC, 1C, A, D, F, VOL, RMT, BC, QF, 

1 TAU, CC ,i CHB, CH, CPH, OQ, V-, VI, CPB, CPT, CADK, X, DT, 

2 IRK, LOC. ITERX, INDEX, I, J: M, N 

C nXYHFwor.LCBIN - PH EMPIRICAL FUNCTIU'I 

C EQUATION 3, A , 

PCF2 «tU0.00668l5*ZZ>~0.1009B)m » 0.449211**2-0. ASA 
RETURN 
■ ENO 

FUNCTION RCr3(KK) 

DIMENSION CIAOI, XNIA0.2), SV!18,50), VTRANU 6 } , RK (1 A, A) , 

1 SCI 14,5! , DCUA), A(6), 0115), F(20), VOLIIO), RMTC2) » 

2 BCIA), QFI6), TAUI5), CC(3>, CHB 1 3) , CH(A), CPHI3), 

3 DO I A) 


1 

03 
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COMUQS/2/ C, XN, SV, VTRAN, RK, SC, DC, A, D* F, VOL, RMT , BC, OF, 

1 TAU, CC, CHH , CH, CPH, DQ, VE, VI, CP8, CPT, CADK, X, OT, 

2 IRK , LOC, 1 T FRX , INDEX, I, J, M, N 
VTRAN FUNCTION 

VARIABLES WITH TIHE DELAYS USED IN EQUATIIHS 8.1-8. 1 . 

RCF3 = SVIKK.LOCI ♦ (SVIKK.LOC * 1) - S / 1 KK,L0C) >PDT/D< 14) 

RF TURN 
END 

FUNCTION SSC2WIX> 

C CALCULATION OF ‘STEADY-STATE OXYGEN REQUIREMENTS FOR VARIOUS LEVELS 
C CF WORK LCAO IX=V)ATTS). 

SS02W = .195 ♦ IX/84.15) 

IF ( X .GT. 210. J SS02H => 2.7 

RETURN 

END 

FUNCTION SS VENT I X I 

C CALCULATION OF STEADY-STATE VENTILATION RA 'E AS A FUNCTION OF TISSUE 
C OXYGEN METABOLIC RATE. 

ir(X.LE..195) SSVENT =5.398 
IFIX.GE.2.)SSVENT=5S. 364-50.*lX-2.) 

IFI IX.GT..1951 .AND.I X.LT.2. ) JSSVENT«27.u3*X 

RETURN 

END 
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